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DESIGN PROCEDURES FOR DOMINANT TYPE SYSTEMS 
WITH LARGE PARAMETER VARIATIONS 

Abstract — This report presents design procedures for 
fourth-order dominant type systems with large plant 
parameter variations. The s-domain specifications of 
the system are assumed to be in the form of an accept- 
able dominant closed loop pole region and bounds on 
the location of the "far-off" closed loop poles. 

The design philosophy is to place compensation zeros 
within the acceptable dominant closed loop pole region 
such that the dominant closed loop poles remain within 
their prescribed region despite the large variations 
in the plant parsmeters. Design procedures are pre- 
sented for variation in the plant gain factory only 
and for simultaneous variation in the plant gain 
factor and the plant poles. Finally, an approxi- 
mate procedure is presented which considers simul- 
taneous variation in the plant gain factor, the 
plant poles and a plant ^ero located on the real 
axis in the s-plane. In all cases, the design 
procedures are such as to minimize the sensitivity 
of the system to internal noise. 
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CHAPTER I 

PROBLEM STATEMENT AND DESIGN PHILOSOPHY 

1 . 1 Statement of the Problem 

The problem considered in this document may be 
stated as follows: l). A single input-single output 

plant* P(s) has parameters (gain factor, poles arid 
zeroes) which may lie or "slowly” vary within a given 
region in the s-plane. 2). The acceptable region of 
closed loop poles is specified in the s-plane. This 
acceptable region is determined by transforming the 
typical type domain specifications such as rise time, 
over-shoot and settling time into specifications on 
the location of the dominant poles of the system. 

This transformation of these specifications fi’om the 
time domain to the pole-zero domain is considered by 
Barber. ^ 3^- Linear time invariant compensation is 

to be chosen to satisfy the above specifications such 
that the effect of internal noise at the plant input 

* 

The term "plant" is used here in the commonly 
accepted sense in the control literature to denote the 
constrained part of the system whose output is ohe 
system output . 



2 


is minimized. The pictorial representation of a typi- 
cal problem encountered in flight control is shown in 
Fig. 1.1, 

1.2 State of the Art . „ 

2 3 

This problem has been treated in the literature. 
However, the design techniques presented have several 
deficiencies. The mapping of the region of plant 
parameter variation into' the acceptable closed loop 
pole region is approximate and is only valid if the 
acceptable closed loop pole region is relatively small 
-? and well removed from the plant parameter variation in 
the s-plane. An additional deficiency is that the un- 
avoidable large loop transmission bandwidth involved 
in using these design techniques results in a very 
unfavorable high frequency response to internal noise 
at the plant input. 

A design procedure which greatly alleviates 

these deficiencies is the topic of a recent paper by 

4 

Horowitz. The design procedure presented in the 
present document is essentially that developed in the 
paper by Horowitz, expanded to take into account the 
effect of one of the "far-off" poles. Thus Horowitz 
uses a third order representation in deriving the 
dominant poles of the system whereas this treatment 
uses a fourth order representation. Chapter IV.ilso 
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considers the additional effect of a drifting zero 
on the real axis. , 

1 .3 Design Philosophy 

The design procedure presented in this paper 
is based on the dominant pole concept. Compensation 
zeroes are strategically placed near or within the 
acceptable closed loop pole region (See Figure 1.1). 
The gain of the system is then determined such that 
the closed loop poles remain within the acceptable 
region in the s->plane , despite the large plant para- 
meter variations. "Far-off" poles (and zeroes) are 
then assigned in such a manner as to allow the loop 
transmission to decrease as fast as possible without 
violating the system time domain Specifications. The 
problem of the placement of these "far-off" poles and 
zeroes, excluding the nearest "far-off" pole placed on 
the real axis, is not considered in this paper. 

Zj. 

Horowitz has presented a method for the placement of 
these "far-off" poles and zeroes. 

1.4 Scope of Work and Terminology 

The plant parameter variation and acceptable 
dominant closed loop pole region are somewhat similar 
to those encountered in flight control. 
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Chapter II considers the problem of variation 
in the plr.nt gain only. Plant pole cancellation and 
replacement is anticipated. The plant transfer 
function P(s) is assumed to be of the form 


P(s) 


k 

2 

s(s +S s+P ) 
v P P' 


( 1 . 1 ) 


where: k = gain factor of the plant which may vary 

between k . ' and k ; 

mm max’ 

S and P are fixed parameters that determine 
P P 

the position of the plant poles. 

For complex plant poles located at a ± jui , & and P ■ 

P P P P 

are given by 

( 1 . 2 ) 

(1.3) 


S = -2a 
P P 


2 2 

P = a + w 
P P P 


The loop transmission L. (s) is assumed to have the fo 


rm 


kk(s +S s+P ) 
«--(•) ■ — - ° - ° 


khn,(s) 
u d ' 1 


S ( 8 +S^ S + P^ )(s+P 1 ) 




(1.4) 


where: K = fixed gain added to the system; 

S q and P & are fixed parameters that determine 
the position of the compensation zeroes; 


.1 ,*• ; r ■ 
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S. and P fl are fixed parameters that determine 
Jf l * 

the position of the plant replacement poles; 
-P^ is the location of the nearest "far-off" 
open loop pole on the real axis. 

For complex compensation zeroes located at a ± juu , S 

Z Z u 

and P q are given by 


S = -2a 
o z 


(1.5) 


„ 2 2 

P = a +uj 
o z z 


( 1 . 6 ) 


For plant replacement poles located on the real axis 


at r^ and r^ , and are given by 


S l = -( r l +r 2> 


(1.7) 


P f ~ r 1 r 2 


( 1 . 8 ) 


The closed-loop transfer function T^(s) is assumed to 
have the form 


P r P f P f P r P f P f 

T / v 1 2 A 1 2 

® ) ~ 2 — n ( o ) 

(s +S s+P )(s+p )(s+p f ) d ' ' 

r r i 1 i 2 


(1-9) 


where: -p*. , ”P#« are the positions of the non- 

' ... 1 . - *2 <'■ : : : / 

dominant closed loop poles which \ 

'T »ay be real or complex conjugate; ' 

V 3 r and P r are parameters that determine the 


F o r c c omplex 


of the dominant closed loop poles . 
closed^! loop polesv^oca^ed^^t- 


rv_ 

a 


JP .. : are given by 


S r > ^o d 


2 2 

? r f 


(tiiaX 

w r ' " " rr _?• 7- ■" ■ 

; (l.tli 


- Chapter XIX considers the problem of simultane- 
ous plant gain and plant pole variation. The plant 
transfer function has the same form as fn Eq . I .=1 but 
S_ andP_ are now slowly varying or unknown parameters 
Since plant pole cancellation is not practical here-, 
the- loop, transmission has the form c " 


ki.(s 2 +S fP ) 

t / \ : O; f » / 

L d<-> - T5 


S(S t S- rP p ){s+P,) 


4 kXn d (s) „ ... 

- d„(») [' - ,2) 


The : closed loop transfer function is. of the same form 
as in Eq. 1.9/ ' . . c 

Chapter IV considers the -problem of simultane- 
ous plant gain, plant pol-e-and plant zero variation. 
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The plant transfer function is assumed to be of the 
form 


l>(s) 


k( s+z ) 

s(s 2 +S s+P ) 
P P ' 


( 1 . 13 ) 


where: -z is the position of the drifting zero on 

the real axis . 

The loop transmission has the form 


L d (s) 


klv( s^+S^s+P^) ( s+z ) kKn^(s) 

s(s^+S s+P )(s+P ) d d ^ S ^ 

v p p' v z 7 


( 1 . 14 ) 


where: -P is the position of the fixed pole used to 

partially cancel the effect of the driftirg 
zero. 

The closed loop transfer function is 


T d (s) = 


P r P f P A (s.z) 

1 z _ 

2 

(s +S r s+P r ,) (s+p f )(s+P c ) 

1 z 


where : 


P r P f p c //z ( s+z ) 

1 z 


-p is the position of the 
c 

z 

near the drifting zero z on 


-p is the position of the 
1 

loop pole on the real axis. 


( 1 . 15 ) 

closed 1 . j p pole 
the rea 1 ax> s ; 


far-off " c : osed 
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Chapters II and III are extensions of work by 

k 

Horowitz in that the effect of the "far-off” pole 
has been included while Chapter IV is completely 
original . 

Appendix A presents various convergence proce- 
dures for factoring polynomials on a digital computor. 
Appendix B presents a geometric proof of a design 
procedure used in Chapter II. 




Cli AFTER II 

DESIGN FOR PLANT GAIN VARIATION ONLY 

2 . 1 Plant Pole Cancellation 

This chapter treats the case when the plant 
pole variation is " sufficiently email " such that 
cancellation of the plant poles is valid providing, 
of course, that they are not located in the right half 
plane. A discussion on what constitutes "sufficiently 
small" variation in the plant poles is given in (4 j. 

The saving in the gain-bandwidth product of the loop 
transmission, obtained by cancelling the plant poles 
and replacing them with poles nearer the desired 
closed loop pole region, may be quite substantial , 
especially if the acceptable closed loop pole region 
is far removed from the original plant poles. 

2 . 2 Design Phil o sophy 

The design philosophy in the case of plant gain 

variation only is to first cancel the existing plant 

poles and replace them with poles nearer the desired 

closed loop pole region. These poles will be near or 

on the boundary of the acceptable closed loop pole 

region at minimum plant gain, k = k . „ 

min 


Compensati on 



zeroes are then located so that the dominant closed 


1 1 


loop poles lie within the acceptable region despite 
the variat ' 'ns in plant gain factor k. The problem 
is depicted pictorially in Fig. 2.1. 


2 . 3 Design Equations 
The expr 
is from Eq. 1 . 4 


The expression for the loop transmission L^(s) 


kk(s +3 s+P ) 
* / \ V O o' 

L,(s) = 


A ^ n d (s) 

d s ( s 2 + S^ s+P^ ) ( s + P 1 ) d d ^ 


( 2.1 ) 


and the expression for the system transmission ( s ) 

is from Eq. 1.9 


T d (s) 


P r P f P f 
1 2 


P r p f p f 
1 2 


(s 2 + S r s+P r ) (s + p f ) (sn-p f ) 

1 2 


D H ( S ) 


( 2 . 2 ) 


The characteristic equation of the system is then 


D h( s ) = d H( s ) +kKn n( s ) 


(2-3) 


Equating the zero degree coefficients in Eq. 2.3 gives 


kiu’ = P f P f P r 
1 *2 


o 


(2.4) 
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Let P*. p* 

1 

Also define 


and p* denote P , p_ and at k 

-2 1 1 -2 


k . . 
min 


k 



k . k 
min 


P*p* 
r f, 



(2.5) 


o 

It is very desirable to minimize the value of , which 
is the system gain necessary to bring the root loci to 
the acceptable region of the dominant closed loop poles. 
Also, the high frequency asymptote of.L^s) (Eq. 2.1) 
is hj/s , which is an important factor in determining 
the effect of internal noise at the plant input. ^ Large 
increases the possibility of plant saturation by in- 
ternal high frequency noise. The next section deals 
with the choice of the parameters in Eq. 2.5 in which 
the minimization of k^ is the prime objective. 


2.4 Choice of Design Parameters to Minimize System 
Gain 

The choice of the position of the dominant closed 

loop poles at k = k . depends somewhat on the shape of 

the acceptable region for the closed loop poles. Let 

the position of the dominant closed loop poles at 

k = k . be denoted by 
min 7 

Pd* p£ = a § * ( 2 * 6 ) 


Note that 


1 4 

P r * l p d! 2 = ( a S) 2 + (‘ JU S) 2 ( see Fi e« 2,2) (2.7) 

It is desirable to minimize j p* | and thereby P* f 
since, from Eq. 2.5 this tends to minimize . This 
implies that p* should be located on or very close to 
the boundary of the acceptable region for the dominant 
closed loop poles. 

The next problem is to choose the values for 

p* and p* , the non-dominant closed loop poles at 

k ss k . . These closed loop poles will lie on the 

mm 

real axis for small values of gain since is assumed 

to be on the real axis. The values of p* and p* 

*1 r 2 

will have to be determined from considerations of the 

time domain specifications of the system transmission. 

They should be chosen as close in as possible, since 

from Eq. 2.5, this v»ll tend to minimize but if 

they are too close to the origin, the system response 

can no longer be characterized by the dominant pole 

pair. The latter consideration determines the minimum 

distance (from the origin) of these poles. Henceforth, 

it is assumed that p* and pj. are known. 

*1 1 2 

The position of the compensation zeroes Z and 
Z is now considered. Denote the position of the 
zeroes Z, Z as 

Z,Z = O z ±jU' z 


( 2 . 8 ) 
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Note that 


1 6 


P o = |Z | 2 = a z +(JU z ( Sec Fi S- 2.2) (2-9) 

Therefore jz| should be made as large as possible to 
maximize P q which, from Eq. 2.5, will tend to minimize 
K.^ . If the variation in plant gain is very large, the 
dominant closed loop poles will, at k s k max , be very 
close to the compensation zeroes. This means that 
these zeroes must be located near the boundary of the 
acceptable closed loop pole region. 

2 . 5 Positioning of Dominant Closed Loop Poles and 
Compensation Zeroes 

The method used in this paper to fix the 
position of the compensation zeroes is the same as 
in reference (4). This method is to demand that the 
angle of departure of the root locus from the dominant 
pole p* for k > k . be within a prescribed sector, 
given in Fig. 2.3 as H^p*!!^. The choice of this 
sector is somewhat arbitrary but should be fairly 
general and easily applied to different acceptable 
regions of dominant closed loop poles . 

In Fig. 2.3, the angle of departure, 0 d> of 
the root locus from p* is given by 


$ d = i8o° - (* 1+ $ 2 + 9 o^-e z ) 


( 2 . 10 ) 


where; 4*^ = L - p£ p*j 

*2 • L - Pf/5 

»°° * L p d P d 

e, . l z P . 

*2 " L 2 p 3 

e z . e, * e 2 

Solving Eq. 2.10 for 0^ gives 

e z = ^) d +(^ 1 +$ 2 + 90°)-180° (2.11 


Denote the extreme values of the departure angle 
$ d , which lie within the sector H^p^H^, as 

(/) . . £ fl) . £ j/j . (2.1 

drain w d v dmax N 

The extreme values of 0^ are then 


Zinin 

O 

O 

ON 

II 

♦ V^dmin- 

O 

O 

00 

( 2.13 

Zraax 

= 90° 

* V*V*dma* 

1 

L 

00 

0 

0 

(2.14 


The locus of zero posit .ons, Z, Z, such that 0 Z i 


constant, is an arc of a circle drawn through the 
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points p*, and a third point X on the real axis 


defined by the equation 


Xp d = 5 z /2 


(2.15) 


The proof of this statement is given in Appendix B. 

The design procedure is, then, to locate two 
points and X 2 on the real axis corresponding to 


8_ . and g. 


in Eq. 2 .15* Circular arcs c n and c. 


Zmiri w Zsoax — — ----- ~ ~2 

are drawn through the points p*X-p* and ptX^p*, as 

shown in Fig. 2.3* Locating the compensation reroes 


between the arcs c. and c 2 trill then insure that the 


angle of departure of the root locus from the dominant 
pole p* is within the specified sector H ^ P^p * 


When k is much greater than k , , the closed 

max rain* 


loop poles will, at k ~ be very close to the 


compensation zeroes. Therefore in order to insure 
that the root locus remains in the acceptable region 


as k approaches ^- max ? the angle of entry of the locus 


into the complex zero should be checked. The angle of 
entry £ ^ is given by 


$ e - a i *0-2*^ 3 +fx L -90°-180 c 


( 2 . 16 ) 


where: =- £ -p£ Z 


f ' 


*2 - ^ -pL z 
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a 3 ~ £? d Z 

°k = ^d Z 
90° > LZ Z 


These angles are shown in Fig. 2.4. 

The value may be computed from Eq. 2.5 once 

the zeroes have been located. It should be noted that 

if k is not much greater than k . , the zeroes may 

max mm* 7 

not have to be located within the acceptable closed 
loop pole region. 


2 . 6 Open Loop Poles of L d (s) 

The last step in the design procedure is to lo- 
cate the open loop poles of L^(s). From Eq. 2.3, the 
expression for D*(s) is 

D*{s) = d £j (s)+iv 1 n d (s) 

Therefore 

d d ( s ) = ( s ) -X, n d ( s ) ( 2 . 17 ) 

or 

d d ( s ) - (s 2 + S*s+Pj;)(s+p* )(s+p* J-Xj (s 2 + S q s+P o } 

1 2 


( 2 . 18 ) 
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£.4 CALCULATION Of THE ANGLE Of ENTRY 






=s(s J +(p*^+p*^+S*)s + (p^P|> 2 + s ^(p|- i + P^ 2 )+ p J- k 1 )s 

+ (S*p^p* 2+ P*(p^+p| 2 )-K 1 S o )} (2.19) 

= s(s 2 +S l s+P l )(s+P 1 ) (2.20) 

2 

The roots of (s +S^ s+P^ ) may be complex or real de- 
pending on the location of the acceptable closed loop 
pole region and the value of . The anticipated root 
locus, for a design of this type, for the given accept 
able closed loop pale region, is shown in Fig. 2.5. 


2 . 7 Design Example 

The following is a 4 * application of the design 
procedure for gain variation only. 

The de s ign equa t i on s a re 


T §( S ) 


P* 


( 8 +S*s+P* ) ( s+p* )(s + p* ) 
*1 z 2 


k . K(s 2 +S s+P ) 
ram v o o 7 

s(s 2 +, c ^ S+P^ ) ( S +P-| ) 

where T^(s) and L*(s) are the system transmission and 

loop transmission respectively, at k = k . . The 

J ’ min 


L M S ) = 


plant transfer function is assumed to be of the form 


s(s +S s+P ) 
v P P 


where k y the plant gain factor, may vary from 
to k max =1000. The complex plant poles are assumed to 
be fixed and have been cancelled by zeroes placed 
near them. 

The closed loop pole values -p* and -p* are 

1 1 Z 2 

assumed to be -10 and -15 respectively. These choices 

for p|^ and pf^ are purely arbitrary other than that 

they must lie to the left of the boundary shown in 

Fig. 2.6. This report does not show how these values 

of pf* and p£ are obtained but how the design pro- 
1 2 

cedes once they are known. The acceptable region for 
the dominant closed loop poles, as well as other per- 
tinent quantities for this design example, are shown 
in Fig. 2.6. 

As a first choice for p*,the dominant closed 

d 

loop pole at k = k . , let 

* min’ 


P* = -3+ j3 

since this is close to the minimum value of P for this 

r 

acceptable dominant closed loop pole region. The angle 
of departure of the root locus from the dominant pole 
p* must be within the sector H^pJH^* i.e. 


1 40° s ib, * 230° 
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The angle of departure of the root locus from the 
dominant pole p* is given by 


*d - 

18O°-(^ 1 +^ 2 +9O°-0 z ) 

where: = 

i-Pf P* d - Tan*' | = 25.<*° 

*2 * 

L - Pf 2 PS - Tan ' 1 T2 “ 

8 Z = 

^zp*t Lz?* 


From Eqs . 2.13 and 2.14, 


8 - . are given by 
Zmin 


the values of 02max anc * 


e„ = 6 j _qo° + 4> +4> 

Z max v dmax 1 2 

= 230 o -90°+25.4 o +14.0° 

= 179.4° 

fl . = »h . -90°+$,+v 

Zmin • dmm 1 2 

= i40 o -90°+25.4°+i4.0° 

= 89.4° 


Next, two points on the real axis , X 2 are deter- 
mined such that 


Lx i p d = Z f a -~ = 90° 
/_X 2 P* = 9z ”- n = 45° 


2 
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Two arcs are now drawn through the points p*X^ p* and 
P*X^P* (See Fig. 2.6). The construction begins by 
locating a point on the real axis that is equidistant 
from the points p*, X and p*. Using this point on the 
real axis as the center of a circle, circular arcs are 
drawn through the points p^, X and p* with a compass. 

As a first choice for the position of the 
compensation zeroes Z arid Z, let 

Z, Z = -6 ± jO.5 

This position for the compensation zeroes will tend to 
maximize P for the circular arcs c. and c 0 shown in 
Fig. 2.6. Large from Eq. 2.5 will mean a smaller 
value for X^ which is the object of this design 
procedure. 

The expression for n^(s) is 

n d (s) = (s-z)(s+z) 

= (s+6* j0.5) (s+6- j0.5) 

= s 2 +12s+36.25 


The value of X^ may be computed from Eq. 2.5, i.e 


X = k . X = 1 *X 
1 trim 

*1 *2 


_ (.1.8 1(10) Q?) 

36.25 


= 74.5 
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From Eqs. 2.17 and 2.19. the expression for d^s) is 

d d ( s ) = ^d^ S ^ ~^1 n d^ s ^ 
or 

d d (s) = s{s 3 +(p*^p*^ + Sj)s 2 +(p*^p*^ + S*(p|.^+p|.^)+P r -h 1 )s 
+ ( s * P J iP * 2+ Pr(pJ i+ P| 2 )-Ki S o)} 

= s{s 3 ( 10+1 5+6)s 2 +(l 50 + 6 ( 25 )+1 8 - 74 . 5) S 
+(6(150)+18(25)-12(74.5))} 

= s(s 3 +3l s 2 +243-5s+456) 

To obtain the open loop poles of L^(s), this equation 
must be factored. Using the methods in Appendix A, 
this equation is factored into the following open 
loop poles 

d d (s) = s ( s 2 +28 . 4s+ 1 68) ( s + 2 . 65 ) 

= s ( s+ 2 . 65 ) (s+ 9 . 00 ) ( s+ 1 9 • 44 ) 

f rum which 

P 1 = 19.44 

(s 2 + S^s+P^)= (s+2 . 65 ) (s + 9. 00) 

= s 2 +1 1 .65S+23.85 

The angle of departure of the root locus from the 
d ominant pole p* is 
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0 d 180°-($ 1 +$ 2 +90 0 -9 z ) 

0 Z = Z.Zp* + <LZp* 

=■ Tan -1 + Tan -1 2^2 _ 39 . 80 + 49 . 3° 

= 89.1° 

Therefore 0 d = l80°-(25 ,4°+l4 . 0°+90°-89. 1 ° ) 

= 139. 7° 

which is satisfactory since it is very close to the 
minimum value of l40°. The angle of entry of the root 
locus into the zeroes is given by Eq. 2.15, i.e. 

= a 1 +a 2 +a^+a^-90°-i8o° 

= Z--P* z+ /L-p* Z + Lp$z+ Z.p*Z- 270 ° 

=, Tan -1 -^l+Tan' 1 • i -|+Tan" 1 _ ^§-+Tan _1 ^^-270° 
= 7 . 1 °+3.2°+2l9. 8°+130.7°-270° 


= 90 . 8 ° 


which is 
For this 
of entry- 
factory * 
is shown 


satisfactory for this design (see Fig. 2.6). 
particular placement of the zeroes, an angle 
between 20° and 160° would probably be satis- 
The approximate root locus for this design 
in Fig. 2.7. 


2 . 8 Improvement in Design Example 

The first design could be improved. If the 
compensation zeroes could be moved further to the 
left, it would, from Eq. 2*5, decrease the value of Kj . 
With this as the objective, let the second choice for 
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P*l be p* = -3 . 5+ j2 .0 


The angle of departure is constrained to be within 
the sector defined by the equation 

110 ° ^ ^ 230 ° 

Quantities pertinent to this design are shown in 
Fig. 2.8. 

The angles and are 


*1 

= Tan -1 

2 

6.5 

= Tan -1 0.308 

= 17.2 

*2 

= Tan"^ 

2 

11.5 

= Tan -1 0.1 74 

= 9.85 


The same method as before will be used to fix the 
position of the compensation zeroes. The angle of 
departure of the root locus from p* is 

0 d = 18O°-(«' 1 +4> 2 + 9 O o - 0 z ) 

= 90°-17.2 o -9.85 o + 6 z ) 


or 

Therefore 


0 = 0-62.95* 


0 Zmax = 230°-63° 


0 „ . = 1 1 0 °- 63 ° 
Zmxn 


n 

47° 


Again two points on the real axis X^,X^ are determined 
such that 

£X lP . = 83.5° 


ix 2 pj = 23.5" 





Two circular arcs are now 


drawn through the points 

p* X* p* and p*X«p* as shown in Fig* £ . § . 

*g l*d *a 2*d 

The compensation zero posit ions a re chosen as 
l t Z - -7± j t 


The expression for n^(s) is then 

P 

n d'~^ K s *+l4s+50 
The value of from Eq. 2.5 is 

P rPF P f 


K 1 " 


1 2 _ ( 1 6 « 2 3 ) { 1 0 ) •' 15 ) 


50 


= 48.69 


This is a reduction of about 4 db« from the X 1 of 
the previous design. 

The angle of departure £1 ^ of the root locus 

from p*, for this position of compensation zeroes, is 

1 20°. The angle of entry $ of the root locus into the 

e 

ompensation zero 7 is 84 c . Both of these values are 
satisfactory for the given acceptable closed loop pole 
region (see Fig c 2.8). The reason that the angle of 
entry should be checked is as follows: If the compen- 

sation zeroes were placed in the extreme left hand 
corner of the acceptable closed loop pole region, an 
angle of entry greater than approximately 80° would 
be unsatisfactory. An angle of entry larger than 80° 
would probably indicate that the root locus would be 
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outside the acceptable closed loop pole region for 

some value of k between k = k . and k = k . 

min max 

The expression for d d (s) is found froir. Kq . 2.19 

d d (s) = s{s^ + (lO+t 5 + 7 )s 2 + (l 30+7(25 ) + 1 6.23-48 . 69 ) s 
+(7(1 50 ) + 1 6 . 23 ( 25 )- 1 4 (48. 69 ))} 

= s{s 3 +32s 2 +292. 5s+773} 

Factoring this equation results in the following ex- 
pression for d^(s) 

d d (s) = (s+4.46)(s 2 +27.7s+173) 

= (s+4.46}(s+9.49)(s+i8.21 ) 

The root locus for this design is of* the same form as 
the previous design and is shown in Fig, 2.9* 

2 . 9 Summary of Design Procedure 

The design procedure for variation in plant gain 
factor only is summarized below* 

1. Cancel the plant poles. 

2. Fix the position of the dominant closed loop 

poles at k - k . . 

min 

3. Determine the values of the "far off" closed 

loop poles at k = k . . 

4. Determine the position of the compensation 

zeroes so that the root locus from k x; k . 

mxo 
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to k=k remains within the acceptable 
max 

closed loop pole region. 

5. Solve for the open loop poles of L^(s) 
from Eq. 2 . 20 . 

The next chapter in this paper considers var- 
iation in both the plant poles and plant gain factor. 






CHAPTER III 


PROBLEM OF SIMULTANEOUS PLANT GAIN 
AND PLANT POLE VARIATION 


3 . 1 Problem Definition 

In the case where the plant poles vary- as well 
as the plant gain factor, plant pole cancellation is 
not feasible. In this case the system gain must be 
sufficiently high so that the dominant closed loop 
poles remain within their acceptable region despite 
the variations in the plant poles. The problem re- 
solves into locating the compensation zeroes such that 
the system gain necessary to accomplish this is minimized. 

3.2 Design Equations 

The expressions for the dominant part of the 
loop transmission and system transmission given- by 
Eqs. 1.12 and 1.9 are repeated below. 

kK(s 2 +S o s+P o ) kXn d (s) 

s(s 2 + S p s + P p )(s + P 1 ) 




P f P fA 


(s +S s+P )(shp )(s+p ) 

1 *2 


P f p r 
1 2 

“ D d< 8) 


Tim plant, transfer function is of the form 
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P(9 ) _h 

s(s +S s+P ) 

P P 

The plant gain factor k varies from k = k . to 

rain 

k = k raav and the. plant poles may lie or "slowly vary” 
within the region shown in Fig. 3.1. The character- 
istic equation of the system B^(s) from Eq. 2.3 is 


D d ( s ) = 

or (s 2 +S r s+P r )(s+p f ^ )(s+p f ^) = s(s 2 +S p s+P p )(s+P 1 ) 

+kh(s 2 +S o s+P o ) 

(3.1) 

s<, * (S r*' > f 1 ** , f 2 > 53 *<« , f ) Pf 2 * S r(f’f ) *f’f 2 )* p r> s 2 *-< S rPf ] Pf 2 

* P r(< , f | tp r 2 >)^ P r p f t P f 2 * s' , .( s p .P 1 )» 3 t ( S p P ) »P pt kk , s 2 

: +(P P,+kKS ) s+kXP 

v p 1 o 7 o 

(3.2) 

Equating the coefficients of Eq. 3.2 yields the follow- 
ing set of equations 

S r + Pf 1 + Pf 2 * S p +P 1 < 3 * 3 ) 

Pf/f/MPf^Pf^r = S p P 1 +P p +kK (3 * 4) 






M 



*• 


r. 


f 10.3.1 PROBLEM Of SIMULTANEOUS PLANT OAlN 
AND PLANT POLE VARIATION 





S r p f 1 f V P r<' , f/' > f 2 > ■ P p P .* kkS „ 

( 3 . 5 ) 

P r p f 1 P f 2 * 

( 3 . 6 ) 

The following substitutions in Eqs. 3.4 

and 3-5 are 

made 


p f, +p f 2 - V p i- S - 

( 3 . 7 ) 

kivp 

o 

Pf Pf - p 

X 1 x 2 r 

(3.S) 

kXP 

P r ° * S r ( V P 1 - S r> +P r = S p P 1 ♦V' kK 

(3.9) 

S kKP 

% ° +P r(V P '’ S r ) * P p P '*“ 

r 

(3.10) 

Define the following quant ites 


Y = kh.P 
' o 

(3.11) 

X = S P 1+ P +kX 
P 1 P 

(3.12) 

Y = P P.+kKS 
pi o 

(3.13) 

Equations 3.9 and 3.10 are then 


+s (S +P.-S )+P = X 

P r' p 1 r 7 r 

r 

(-.14) 

T*- * P r<V P l- S r) * Y 

r 

(3.15) 


Considering the equations for X and Y, (Eqs. 
3.12, 3.13), the variation in X, &X, and the variation 
in Y, AY, due to the variation in plant poles only, 

i.e., parameters S and can be expressed as 


AX = P- AS +AP 

1 p p 

(3.16) 

AY ^ P ? AP p 

(3.17) 


The variation in plant gain factor will be considered 
later in the design. 

3 - 3 D esign Procedure 

An outline of the design procedure that will be 
followed > n this problem is as follows: : 

1. ) Map the acceptable region for the don.^nant 

closed loop poles into the X,Y plane using 
Eqs. 3.1^ and 3.15 for fixed values of the 
parameters y * and S • 

2. ) Map the plant pole variation into the AX, 

AY plane using Eqs. 3.16 and 3.17 for fixed 
values of the parameter P^ . 

3. ) Compare the two mappings in (1,2) above. 

If the mapping of the plant pole variation 
in the AX, AY plane does not fit ir to the 
interior of the mapping of the acceptable 
dominant closed loop pole region, the map- 
ping of the latter will have to be repeated. 
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using a larger value of y . 

4.) Solve for the values of kK, S and P 
using Eqs. 3.11, 3*12 and 3.13, where 
the values of X and Y are obtained from 
the positioning of the mapping of the 
plant pole variation in the interior oi 
the mapping of the dominant closed loop 
pole region in the X, Y plane. 

The next four sections in this chapter elaborate on 
these four steps in the design procedure. 

3.4 Mapping of the Dominant Closed Loop Pole Region 
The mapping of the acceptable dominant closed 
loop pole region into the X,Y plane involves the para- 
meters y, Pj and Large y implies large gain, i.e., 

large kk, since P q does not have a large range of 
values (See. Eq. 3.11)* An approximation for the 
value of y can be obtained from Eq. 3.6, i.e. 


Y 


Pf Pf i 
I 1 I 2 


(3.18) 


The maximum value of P^ can be found from the acceptable 
r gion for the dominant closed loop poles. Denoting 
the value of the dominant closed Loop pole by = 
a^+ju>d> the value of P^ is 


P 


r 



a 



(3.19) 


39 


"'ho values of p arid p can be roughly approximated 

1 1 2 

by considering the boundary for these closed loop poles 

shown in Fig. 3*1. These poles will be complex for 

large values of system gain. 

, the nearest h far off " open loop pole located 

on the real axis , should be chosen as close in as 

possible, since from Eq. 3*3 , this will decrease the 

values of p r and p which, from Eq. 3-8, will tend 

to decrease the value of fixed gain that must be added 

to the system. If is chosen too close in, though, 

the closed loop poles p and p„ may lie to the right 

1 1 1 2 

of the vertical boundary shown in Fig. 3.^ violating 

the specifications of the problem. 

' The major problem in this mapping operation is 

the parameter S . Denoting the value of the plant pole 
P 

as p = a +juu , the value of S is 
P P P 


S = 
P 



( 3 . 20 ) 


In the point by point mapping of the acceptable region 

for the dominant closed loop poles, there is no criteria 

for determining what value of S to associate with a 

P 

particular point on the boundary of the acceptable 
dominant closed loop pole region. This dilemma is re- 
solved by considering the following argument: If the 

dominant closed loop poles are to lie 'within their 


acceptable region despite the variation in the plant 

poles, the mapping of the dominant closed loop pole 

region can not be highly sensitive to the position 

of the plant poles and hence the value of The 

actual point by point mapping of the dominant closed 

loop pole region into the X,Y plane is performed using 

several different values of S for each value of y anci 

P 

P- • These different values of S should include the 

1 p 

minimum arid maximum values of S for the given region 

P 

of plant pole variation as well as values in between. 

A median value of S is denoted by S , in Fig. 3*2. 

p 7 pmed 

The mapping of the acceptable closed loop pole region 

into the X ,Y plane is depicted graphically in Fig. 3*2. 

The values of p„ and p r , . the " far off" closed 

loop poles, are also of interest in the mapping of the 

dominant closed loop pole region. The values of p and 

i 1 

p~ must lie to the left of the boundary shown in Fig.3*1 
2 


For each point on the boundary of the 
loop pole region and given values of 
y, and S^, the values of p^ and p 
as follows; From Eqs. 3*3 and 3*6 


dominant closed 
the parameters 

may be obtaine 

*2 


d 


% +p e 2 = V p i‘ s r 


P f P f = Y/P r 

1 l 2 


( 3.21 ) 
( 3 . 22 ) 


= S , 


Define 


P 


+ P 


(3.23 ) 


hi G-. 3.1 MAPPING OF THE DOMINANT CLOSED LOOP 
POLE REGION IN THE S- PLANE INTO THE Y 


Plane 


MAPPING- EQUATIONS 
X = X +S*(Sp*P A -$ i ) 


V * Lh. + fit (£* -S A ) 

PARAMETERS 

>t IV 

S- PLANE X-Y PLANE 


R. ~ ■* 


C^AL>< 



u, y) 



3 

s p * s /»/f«r 
5 /* * 


ACCEPTABLE CLOSED 
LOOP POLE PEG-iON 




■V'2 = " f 


living bqs, 3 *23 and 3«2** Tor p 0 arid ;\eldt 

1 d 


= * /(S,/2) J 


f* 2 


% = -r - 


(3.25) 


(3.2o) 


If the values of p r and o r fall to the right of the 

l 1 ‘i 2 

vertical boundary shown in Fig* 3*1* P- has been placed 

too far in. The values of and obtained during 

t 1 1 2 

this mapping operation will not correspond exactly with 
those in the final design, since the boundary of - the 
plant pole variation will not, in general, map exactly 
onto: the boundary of the acceptable dominant closed 

loop pole region. Additional features of this mapping 
operation are covered in Section 3*8* 


3 . 5 Mapping of the Plant Pole Variation 

The only parameter in the mapping of the plant 
pole variation into the £X, L Y plane is P^ (See Eqs. 

3.1o and 3.17J* The value of used in the mapping 
must be the same as that used in the mapping of the 
dominant closed loop pole region -Eqs. 3-1**, 3- 3). The 
mapping of the plant pole variation 3 s it.pi einen - ed by 
defining any point (°po* x> p 0 ) on boundary of "die 

plant pule variation. The nominal values ol b and P 


plant pule variation. 


and P 


4 


are then 


(3.27) 


Vj 


s 

po 


V 

pt> 


-2a 

2 

o 

po 


po 

2 


po 


(3.28) 


Equations 3.16 and 3.1 J may then be written as 

AX = P t (S -S ) + (p -F ) (3.29) 

1 v ppo /x ppo / x ' 

AY = p i (P p-‘ > po^ (3.30) 

The region of plant polo variation is then mapped 
point by point into the ^X t AY plane, as shown in 
Fig. 3*3. It should be noted at this point that the 
shape and size of the mapping of the plant pole varia- 
tion in the aX, AY plane is not dependent on the choice 

of the point (a ) and hence on the values of S 

' po f po' po 

and P . Different choices for this point will only 
alter ti e po rtion of the mapping in the AX, AY plane. 
The units on the AX. AY axes In the plane must he 

the same as those on the X,Y axes in the mapping *f 
the dominant closed loop pole region in the X,Y plane . 
When the mapping of t..e plant pole variation in the 
AX, AY plane is transferred to the X,Y plane, its 
angular position with respect to the AX, A Y axes must 
be preserved, x.e., the mapping of the plant poJe var- 
iation may not be rotated in the \.Y plane. Additional 
features of this mapping operation are covered in 
Section 9 of this chapter. 


FIG. 3.3 MAPPING OF THE PLANT POLE VARIATION 
REGION IN THE S-PLANE INTO THE PLANE 

MAPPING EQUATIONS 

AY. 

parameter. 5 
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The problem now is to fit th^ mapping of the 
pla nt pole variation in the iiX, a Y plane into the 
interior of the mapping of the dominant closed loop 
pole region in the X,Y plane as &hown in Fig. 3.H. If 
this fit is not possible, the mapping of the dominant 
closed -loop pole region will have to be performed for 
larger values of y * 

3 .6 Calculation of the System Gain and the Compensation 
Zero Location 

Once the mapping of the plant pole variation 

fits inside the mapping of the dominant closed loop 

pole region in the X,Y plane, the value of the system 

gain kk and the compensation zero positions, given by 

the parameters S q and may be computed. For a linear, 

time invariant, minimum-phase system, a value of system 

gain kK can always be found such that the mapping of the 

plant pole variation will fit inside the mapping of the 

7 

dominant closed loop pole region in the X,Y plane. 

Figure 3.4 shows the mapping of the plant pole 

variation fitted inside the mapping of the dominant 

closed loop pole region in the X,Y plane. To solve 

for the values of kK, S q and P q , a point on the boundary 

of the mapping of the plant pole variation is chosen 

where the values of S and P are known. In Fig. 

P P 

this point is denoted by A. The values of S and P at 

P P 

point A are, from Fig. 3.3 and Eos. 3.27 and 3.28, 


FlG. 3 .* CALCULATION OF- THE SYSTEM GAIN AND 
THE EERO LOCATION FROM THE MAPPINGS OF THE 
oominant Closed loop pole region and the 

PLANT POLE VARIATION 



r(*spec t ively , S and P „ Denote the coordinates of 
"a p a 

point A in the X,Y plane as X and . Since y and P t 

a a 1 

are known for this particular mapping, kk, S q and P (> 
may be obtained from Eqs. 3*11? 3-12 and 3.13* as 
follows: 


kk = X -S P.-P 

a Pa 1 Pa 


P = 


kk 


t> = 


Y -P P, 
a P a 1 

kk 


(3-31) 

(3.32) 

(3.33) 


The value of KK should be interpreted as the necessary 
value of system gain, i.e. The actual value of 
added gain to the system is, from Eq. 3.31 



(3.3*0 


mm 

The choice of the point used to compute the values of 

KK , S and P has no effect on the values obtained for 
o o 

these quantities, so long as the point is on or within 
the mapping of the plant pole variation in the X, Y 
plane. To prove tnis, a second point B is chosen, as 
shown in Fig. 3.4. The value of system gain using point 
B is, from Eq. 3.31 

kK * v s P *vv (3 ' 35) 

b b 


Now and P are not known but they can be computed 


using Eqs. 3.29 and 3-30, i.e. 


/ 
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>X = X, 


-X = P. (s -S )+(P -P ) 
>a 1 ' p. p ' ' p. p ' 


p b p a p b p a 


aY = Y -Y = P. (P -P ) 
b a 1 V P b P a ; 


Solving Eqs* 3*36 and 3-37 I'o 1, S and P 


(3.36) 

(3.37) 


results in 


Y -Y 

p = ** a +p 

P b P 1 p a 


X b-V< Y b- Y a>/ P , „ 

s = +s 

Pb P 1 1 


Substitution of Eqs. 3.3$ and 3-39 into Eq, 
yields 

v -Y 

kh = X L -X L +X_+ Y b -Y a -S P. ~ 7 c P 

b b a p a i P, p a 


(3.38) 

(3.39) 
3.15 



This equation is identical with equation 3.3^, which 
implies that the value of k X is not dependent on the 
point used to compute it. From Eqs. 3.32 and 3.33, the 
same statement can be seen to hold for S and P . 


3 . 7 Mathematical Explanation of the Design Procedure 

This section presents the mathematical justifi- 
cation for the design procedure presented in Sectiot. 
3.3. Referring to Fig. 3.5, the dominant closed loop 
pole regions in the s-plane are denoted by C and C. 

The mapping function, which maps the dominant closed 
loop pole region into the X,Y plane, is denoted by T. 
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Denote the mapping of the sets C, C for 

S = S C[S ,,S 1 by means of Eqs. 3. 14 and 3.1 5, 

as S i i . e „ : C,C ■* (See Fig. 3o). Now from 

Eqs. 3.1^ and 3.15, it is evident that no two pair of 
points in C,C map into the same point in S^, therefore 
is a one-to-one mapping function. Since every ele- 
ment of S^ appears as the image of at least one pair 
of points in C,C. T. ma .s 0 and C onto S. . Now. sii ; ce 
Tj, is a one-to-one mapping function and also maps. C 
and C onto , then the inverse mapping function ^ 
exists and maps S. o nto C>C in a one-to-one fashion. 

The trapping of C and C into the X,Y plane for 

an infinite number of values of S between S . and 

p pmxn 

will result in the set shown in Fig. 3-6. This 


pmax 


set may be represented as 


s = n s. 

i= i 1 


(3.36) 


Under the assertion that T^ ^ exists arid maps ont o 
C , C in a one-to-one fashion, i.e., 


tTUs.-.C.C for S = S e[S . ,S ] (3-37) 

11 p Pi pmin’ pmax J w ’’ 

it follows that 


T? :.?Ls,-»C,C for all S € [S . ,S ] (3.38) 

i i=1 l p. pmin’ pmax J ' ' 


Referring to Fig. 3 . 7 * this means that if the 
given variation in X and Y, denoted by the set P , is a 


FIG. 3.7 INVERSE MAPPING Op THE SET P 
IN THE X-Y PLANE TO THE S- PLANE 
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subset of the set S = f) S. , then the inverse mapping 

i~1 1 

of P into the s-plane will be contained in the sets 
C,C. Letting the set P represent the mapping of the 
plant pole variation, then the solution of Eqs. 3.31, 

3.32 and 3.33 yield the necessary values of kk , and 
P Q to position the set P within the set S which repre- 
sents the mapping of the dominant closed loop pole 
region. In this manner the difficult problem of choosing 
kk, S q and P q so as to insure that the dominant closed 
loop poles lie within their acceptable region in the 
s-plane, is transformed into the less difficult problem 
of determining kk, S q and P q such that the mapping of 
the plant pole variation in the ^X, a Y plane may be 
fitted inside the mapping of the acceptable dominant 
pole region in the X,Y plane. 

3 . 8 Analytic Aspects of the Mapping of the Dominant 
Closed Loop Pole Region 

The relative complexity of Eqs. 3.14 and 3.^5 
require that the actual mapping of the dominant closed 
loop pole region be done point by point using a digital 
computer. In this section, the mapping rf curves of 
constant P and curves of constant S in the s-plane 
, • •iiTt'o the X,Y plane is examined. Curves of constant P 

r 

in the s-plane are shown in Fig. 3.8. The mapping 
equations are 


FIG. 3.8 CURVES OF CONSTANT 
P A IN THE S -PLANE/' 
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FIG. 3.9 MAPPING OF CURVES OF CONSTANT P^ 

Vf /N THE S-PLANE INTO THE X~ V PLANE 
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r +S r(V P r S r )+P r 


X 



Since P is to be held constant, S_ must be eliminated 
r r 

in these two equations. The parameters y» P.j and S p 
will also be held constant in these two equations 
Solving the second equation for yields 


S r « 


Y-P P.-P S n 
r 1 r p 


(3.39) 


v/ p r - p r 

Substitution of Eq. 3*39 into the first equation gives 


Y-P P,-P S 


Y-P_P, -P S 


^'v-TPT^i ( V P « 


After considerable algebraic manipulation 
be placed in the form of 


= X 

(3.40) 

Eq. 3.40 may 


-(•'/P r -P r ) 2 {X- Y /P r -P r - 


(S +P. ) 2 


= {Y..±(S p+Pl )( Y /P r +P r )} 2 (3.41) 

Equation 3.41 is in the form 

-4a(X-h) = (Y-k) 2 (3-42) 


which is the equation of a parabola opening to the left, 
whose vertex is at (h,k) with a focal length equal to 
a. delating these quantities to Eq. 3*41, the vertex 


of the parabola, represented by Kq. 3.41, has tho follow- 
ing coordinates in the X,Y plane 


X * Y/P r ^P r + 


<V P 1> 


\2 


( 3 . 43 ) 


Y = i(S p+ P 1 )( Y /P r+ P r ) (3.44) 


The focal .length of this .parabola is 

(Y/P r -P r ) 2 

a *— j— — 


( 3 . 45 ) 


The position of a parabola defined by £q. 3.41 is shown 
in the X,Y plane in Fig. 3.9. r It should be noted that 
Eq. 3. 43 represents a parabola whose axis is parallel 
with the X axis in the X,Y plane. This is indicated by 
the absence of any terns of the form CXY in Eq. 3.41 
with C = constant. 

Curves of constant S in the 3-pltine are shown 

_ •- r 

in Fig. 3»10. Since S is to- be held constant , P 

r r 

must be eliminated in the two mapping equations. 
Defining 


* S V P l ' S r 


(3.46) 


and solving the second mapping equation for P r yields 


P = x ± /xl. 

p r 2A * V kA 2 


S r y 

A 


( 3 . 47 ) 


Substitution of Eq. 3.4? into the first mapping equa- 
tion gives 











= x 


(3.48) 


After much algebraic manipulation, this equation may- 
be put in the following form - 

2 S *P, v 2 " - , y(A-S ) 2 

x jht - 2As r x T y(s p »p, )*A 2 s 2 . — r-j— - o 

r t ' r 

. V (3.49) 

^ - V ■ 

This equation now is in the form of 

- -aX 2 +2bXYtcY 2 +2dX+2eY+f » 0 (3.5P) 

with -- : 7 


a = 1 
b = - 




c 



which is of the form of a second degree equation in a 

8 

rotated system of coordinates. The type of second 
degree equation which Eq. 3.50 represents can be: deter 
mined by examining the coefficients a, b and c. Prom 
Eq. 3.50 

L 2 f j / * p 1 \ y 2 i 1 

b -ac - - i S A 

r r 

= 4 S 2 A 2 {fe + P. ) 2 -S-A] 


(3-51 ) 



r nH 

6 

It can be shown that the following relationships exist 
between the coefficients of Eq. 3.50 and the type of 
second degiee equation. 

2 

b - ac < O *4 an ellipse 

b 2 - ac s 0 a parabola 
2 

b - ac > 0 ■+ a hyperbola 


Equation 3.51 cari put in the form 
4sjA 2 {(S r+ A) 2 -S r A} 

or 4s 2 A 2 {S 2 + S r A+A 2 } (3-52) 

Now if A is greater than zero, so is Eq. 3.52, i.e., 

4S 2 A 2 <S 2 + S r A + A 2 ) > O (3.53) 

since S r is always greater than zero (closed loop poles 
in left hand plane^ For A to be greater than zero, 
bp, P 1 and S p must satisfy the following 

P, > V S p (3.5<.)_ 

The parameter is negative if the plant poles lie in 
the right half plane. Therefore Eq. 3.5** can he 
written as 

p 1 >S r + l s p ! for s p <0 (3*55) 

This equation is certainly satisfied for most, feedback 
control systems. Curves of constant in the s-piane 
therefore map as hyperbolas in the X,Y plane under the 


‘V> 


mapping equations 3-^ arid 3*1 5* The position of a 
hyperbola defined by Eq, 3-^9 in the X, Y plane is 
shown in Fig, 3*^1* 

The angle which the major axis of the hyper- 
bola is rotated from toe X axis in the X,Y plane is 
. 8 

given by 


. 8 


Tan" 1 ! ^ C ~-l 


± J (c -a )'* + 4 b 2 1 
2b J 


( 3 . 56 ) 


= Tan 


_ 1 |( 1 /4S r A-1 ) (V4S r A-1 ) 2 + f ( A + S r )/(S r A) j 2 


- ( 


A+S 

i 

S A 

r 


■) 


•} 


(3.57) 


To obtain an insight on the magnitude of the angle d , 
the following values are assigned to the parameters in 
Eq, 3.57 which are typical for the design example in 
Chapter II and the design example that will be con- 
sidered in this chapter 

S =10 S = 2 

r p 

P 1 = 30 A = 22 

Substitution of these quantities into Eq . 3.57 gives 


8 


Tan 


• 1 ^ (¥ 880-1 ) ± V 880-1 ) 2 - t .( 32 / 220 )‘ 


- (- 22 ) 

k 220 ) 


Tan -1 13-61 = 85.66° 


1 


00 


For the range of parameters considered in this paper, 
a good approximation to Eq. 3-57 is 


3 


Tan 


i 2S A 

-1 r 


S + A 
r 


(3.58) 


The sign in Eq. 3*57 may be arbitrarily chosen to make 
3 positive. 

The only ocher two parameters in the mapping 
of the dominant closed loop pole region are the 

"far-off" pole on the real axis, and y defined by Eq. 
3-11. The qualitative effects of these parameters on 
the mapping of the dominant closed loop pole region are 
shown in Figs. 3-13 and 3-14* The dominant closed loop 
pole region used for this investigation is shown in 
Fig. 3-12. From inspection of Figs. 3.13 and 3-14, it 
is seen that increasing the value of y both increases 
the relative size cf the mapping of the dominant closed 
loop pole region and its coordinates in the X f Y plane. 
Increasing the value of merely increases the coordin- 
ates of the mapping in the X,Y plane. Figure 3-13 
illustrates the variation in the "far-off closed loop 

poles, p and p , during the mapping of dominant closed 
1 2 

loop region into the X,Y plane. These "far-off" poles 
are obtained using Eqs. 3*25 and 3*26 and the procedure 
outlined in Section 4 during the mapping operation. 
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3.9 Analytic Aspects of the Happing of the Plant Pole 
Variation " ■ 

‘ As with napping tff the dominant closed loop 
pole region, the mapping of* the plant pole variation 
(Eqs. 3.29, 3. 30 ) is mo st eas ily a ccomplished by 
mappings point by point using a digital computer. In . 
this- section, the -mapping of curves of coustant P p , 
constant Sp and coustant u in the s -plane -into the 
AX, aY plane is investigated. '' : 

The mapping equations are c r 


AX 


■ P,AS .+AP = P, (S -S ) + (? -P ) 
.1 P O- 1 V p P 7 ' n D ' 


P Poe 


AY = P.AP = Pi (P -P ) 

1 *p - 1 v p p / 


where S and P- define an arbitrary: point a + 

- P ° -o ' • ’ - v P o P o 

on cthe boundary of the plant pole variation in the 

s-plane. c c ’ - 

Curves of, constant P in the s-plane are shown 

. ' . " ' P - . ... . ; 

iF. F ig * 3* ^60 3*he mapping equations with P equal to 

. . r ’ ; - - - P 

a constant are ■ ‘ . c 


AX = P.(S -S )+P -P 

1 P P 0 P P o 


(3.58)' 


AY = P^ (Pp~Pp ) - constant 


(3,39) 


The mapping of curves of constant P in the s-plane 

‘ . P 

into the X,Y plane is shown in Fig. 3.17. >• 





Curves of constant in the s-nlane are shown 

P 

in Fig* 3-18. Eliminating in the two mapping 
equations yields 

uY = P t AX-P^(S -S ) (3.60) 

f P Q 

Equation 3*60 is in the form 

AY = m AX+b (3.61 ) 

which is the equation of* a straight line ir. the AX r< - Y 
plane with slope id and AY intercept b where 


n> = P v ; b = -P?(S -S ) 
l I PP' 

o 

The mapping of curves of constant S in the s-plane 

P 

into the £X,aY plane is shown in Fig* 3.19. 

Curves of constant a* are shown in Fig. 3.20. 
A point on the line is defined as j + juu ^ ^ The 

mapping equations are 


AX = P- AS +AP 
1 P P 

— P* \ **2j +20” ) + (o’ +(li - “C *—& ) 

1 P 1 P P 


( 3 . 62 ) 


AY = P,AP 
1 P 

„ , 2 2 2 2 v 

= P- (c -x- ) 

1 x 1 p o J 


(3-63) 


Since oj^ is to be held constant, a must be eliminated 
in Eqs. 3.62 and 3.63. Solving equation 3.63 for a yields 

(3.64) 


a = ± 


/H 

/ p. 


AY 2 2 2 

+ cr +cu -x, 

P P ' 
o *o 


Substitution of Eq. 3.64 into Eq. 3.62 gives 






fl(r. 3. ZO CURVES Of CONSTANT cj IN THt 

*'*•*"* 


S+jJ, 


itl»U 
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4X = -2P, {±/fW 


+U) 


2 2 

-cu, -a 1 

P i p 

*o o 


(3.65) 


Rearranging Eq. 3.65 gives 


(AX) 2 - 2|XAY + iAYlf -kP^o ax + 4AY (a -P 1 )-4P 2 (u« 2 -UJ 2 )=0 


1 


o 

( 3 . 66 ) 

This equation is now in the form of Eq. 3*50, the form 
of a second degree equation in a rotated system of 
coordinates. The form of this second degree equation 
may be determined by examining the first three coeffic- 
ients, a, b, and c of Eq. 3*66. For Eq. 3*66 


a = 1 ; b = -yP 1 ; c = yP 2 

and b 2 -ac - yp 2 ~yp 2 = 0 

According to the conditions previously stated, equation 
3.66 represents a parabola in a rotated system of 
coordinates . 

The angle which the axis of the parabola makes 
with the £X axis in the AX,4Y plane is given by Eq. 3-36 
which is 


0 = Tan' 1 { ( c - - a . ) - ^(^- a l 2 . tA b i _| 

( 1 /P^-i)W( 1 /pf-0 2 +4/P 2 

= Tan { -2 /p.j ^ 

= Tan -1 P, 


(3.67) 
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where the sign on the radical has been chosen to make 
6 positive. 

Unlike the mapping equations developed for the 
dominant closed loop pole region, equation 3.66 is not 
overly complex and a transformation may be made that 
will eliminate the AXAY term in Eq. 3.66. The trans- 
formation equations which will transform Eq. 3.66 into 
the following form 

a' (AX' ) 2 +c' (&Y> ) 2 +2d'AX»+2e'AY' +f = 0 (3.68) 


where the primes indicate quantities referred to the 

g 

rotated system of coordinates are 


a 1 

2 2 
~ acos 0 +2bsin6 cos9 +csin 0 

(3.69) 

c * 

2 2 
= asin 0 -2bsin0 cos0 +csin 6 

(3.70) 

d' 

= dcosG +esin& 

(3.71 ) 

e 1 

= -dsin9 +ecos0 

(3.72) 


The trigonmetric functions may be expressed as follows 


COS0 = 


ij 1 +Tan 6 Vi +P^ 


sin3 = 


Tan0 


1 


V 1 +Tan 2 0 v 1 


(3.73) 


(3.74) 


The new coefficients a*, c T , d ! and e', are 
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a ' = 


(2/p 1 )p 1 (yp;)?! 


i+p: 


i+p: 


i+p 


= 0 


(3.75) 


c ' = 


P* (2/P 1 )P 1 »/P* 


1+P! 


1 +P. 


i+p: 




(3.76) 


^P.,*? 2P, (a -P, ) 

IP » P * 

* r n 


d* = 


-2p: 


yi+p; 


yi + p; Vi+p^ 


e ' = 


2P ?°p 2 ( a p- P 1> 20 p (P?.')-2P, 

f_o *o _ o 


v ; 1+p 


Vi+p 2 \/i+pf 


(3.77) 


(3.78) 


Equation 3.68 with a ' =0 is of the following form 


c' (AY 1 ) 2 +2d'AX'+2e’AY’+f = 0 


(3.79) 


which may be put in the form of 


<AY'.f;) 2 = - |f72_) 


2dJ 

c 


(3.80) 


From Eq. 3.^2, this equation represents a parabola in 
AX' ,AY' system of coordinates with its center at 


AX s - 


2d 


_ + i-iT 

' 2d’ c' 


( 3 . 81 ) 


AY' = -e '/c 


(3.82) 


and focal length a of 

d' 


a = — 


Ac ' 


(3.83) 
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In most feedback control systems, is much 
greater than unity. Therefore the following approxi- 
mations can be used for the coefficients c 1 , d* and e f ; 


c ' = 


p ? +1 




(3.84) 


2P. 


d' = - 


= - 2P, 


VlH-Pi 


( 3 . 85 ) 


2a (i+pp-2P 1 


e ' = 


= 2 (a P^l ) (3.86) 


The approximation to Eq. 3-80 is then 


{AY'+2(a P^l)} 2 = 4P 1 {AX« + P 1 (x 2 -& 2 ) + 


(a P.-l)' 

PI 7 


(3.87) 


The position of the parabolas defined by Eq. 3-87 in 
the AX, A Y plane is shown in Fig. 3-21. From Eq . 3-87 
and Fig. 3.21, it is seen that the AY* coordinate of 
the center of the parabola in the AX', AY* rotated 
coordinate system is not dependent on the parameter . 

The only other parameter in the mapping of the 
plant pole variation is P { . The qualitative effect of 
on tne mapping of the plant pole variation is shown 
in Fig, 3.23- The plant pole variation used in the 


FIG. 3.ZZ PLANT POLE VARIATION REGION 


US BO TO ILLUSTRATE THE EFFECT OF Pj on 
THE MAPPING- INTO THE AX “AY PLANE 





mapping >.3 shown in Fig. 3.22. From tig. 3.23, it is 
ohnorved that the size of the mapping of the plant 


pole variation increases with increasing . Since 
the size of the mapping of the dominant closed loop 
pole region is not highly sensitive to P (see Fig. 
3.14}, it is very important to bring in as far as 
possible to minimize the needed system gain kh, 

3-^0 The- 'Effect of Pln^tt. Gain Variation an the design 
- . - It is possible that the don inant closed loop 

P^lts may lie Vi thin their acceptable region for k=k 

- - : *■ y. - - min 

and vet move outside the acceptable region for k>k . . 

Figure 3*2^ illustrates this pq^s^ihii it y:. This un- 
desirable possibilit y can £>e predicted by - consider ing 
the angle of departure of the root locus from ihe domi- 
nant closed loop pole at k = k . and the ap^le of 
entry of the root locus into the compensation zeroes. 
This is the procedure: used by Horowitz j » his paper. ^ 
Figure 3-25 illustrates ""tha dominant closed 

loop pale ? within their acceptable region at k = k . . 

mm 

lor a given plant pole variation t the dominant closed 
leap foies for k -= > . may lie anywhere oh a boundary 
such as the boundary ABCP shovn in fig, 3.25. It is 
usually sufficient to check the angle of departure of 
the root locus for a few points arotirrJ the boundary at 


k 


k . 
*nxn 
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To outline the procedure, consider dominant 

closed loop poles at B and B Tor k = k . and "far-off" 

* mm 

closed loop poles at -p,. and -p~ for k = k . . De- 

* ^ f 2 mxn 

fine the following angles 


L B B = 

90° 

(3.90) 

L- P r B = 

a i 

(3.91) 

L- P f B = 
2 

C 2 

(3.92) 

L ZB = 

*1 

(3.93) 

L ZB = 

*2 

(3.9*0 


The angle of departure of the root locus £ from the 
dominant closed loop pole located at B is then 

<& d = l80 o -(a 1 +a 2 +90°-«> 1 -v 2 ) (3-95) 

It is relatively easy to ascertain for a particular 
mapping what angles of departure may lead to an un- 
satisfactory design. For example, for a dominant closed 
loop pole located at point B in Fig, 3.25, an angle cf 
departure & ^ of -10° would probably be unsatisfactory. 

If there is still some doubt whether the root 
locus remains within the acceptable region for a parti- 
cular dominant closed loop pole location, the angle of 
entry of the root locus into the compensation zeroes 


can be checked. As before consider dominant closed 


loop poles located at B and B and define the following 
angles 


L Z Z 

= 90° 

(3.96) 

L B Z 

= 9 1 

(3.97) 

L B Z 

= 9 2 

(3.98) 

L -P f Z 
1 

= 6 3 

(3.99) 

L -p Z 
*2 

= 8 4 

(3.100) 


The angle of entry of the root locus & ^ to the complex 
zero Z for dominant closed loop poles located at B and 
B is then 

= a 1 +s 2 +9 3 +e i .-‘9o o -i8o° ( 3 .^ 01 ) 

As stated for the angle of departure criteria, it should 
be relatively easy to determine for a given mapping 
what angles of entry are acceptable. Considering the 
case again for dominant closed loop poles located at 
B and B, if the root locus had an angle of departure 
from B of -10° and an angle of entry to the. compensa- 
tion zero Z of +10°, this would probably confirm that 
the root locus is outside the dominant closed loop 
pole region for some value of k > k 

min 
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If such a situation arises in the design, the 
procedure is to increase the value of added gain K 
until all angles of departure and entry are satis- 
factory. The object, of course, is to obtain a design 
with the least amount of system gain klv. 

3 . T 1 Design Example 

A design example is presented in this section 
based on the previous design procedure. The region of 
plant pole variation, acceptable dominant closed loop 
pole region and boundary for the "far-off' closed loop 
poles are shown in Fig. 3.26. The acceptable dominant 
closed loop pole region is the same as used in the 
design example of Chapter II. The plant pole variation 
is essentially the same as used by Horowitz.** The 
boundary for the "far-off" closed loop poles has been 
arbitrarily chosen as shown in Fig. 3.26. 

Following the design procedure outlined in 
Section 3 of this Chapter, the dominant closed loop 
pole region is mapped into the X,Y plane with parameter 

P y and S . 

1 P 

After several computer runs, values for and 
Y of 30 and 20,000 respectively were used for the first 
design. The value of 30 for was chosen because this 
value places the "far-off" closed loop poles very near 
the boundary shown in Fig. 3.26. Using E^. 3.18 with 


FIG, 3.2E S’ DOMAIN SPECIFICATIONS FOR. THE 



I 
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P = 60 and p_ = p_ = +15“j10, a value of approx- 
r *1 X 2 

imately 19»000 for y is obtained which is quite close 

to the value used in the first design. The value of 

for this approximation was found from the acceptable 

dominant closed loop pole region while the values of 

p^ and were estimated from the knowledge of the 

1 ^ 1 2 

anticipated root locus. The values of 3^ used were 
-o, 2 and 10 which correspond to the minimum, median 
and maximum possible values of for the plant pole 
variation shown in Fig. 3*26. As anticipated the 
mapping of dominant closed loop pole region is not 
highly sensitive to the value of at this large 
value of y. The mapping of the dominant closed loop 
pole region for these parameter values is shown in 
Fig. 3.27. 

The mapping of the plant pole variation for 
various values of including =30 is shown in 
Fig. 3.28. 

It should be emphasized that the only practical 
means of performing these mapping operations is on a 
digital computor. On such a machine, the mapping of 
the dominant closed loop pole region may be easily per- 
formed for many different combinations of y and '■? . 

The same values of P^ are then used in the mapping of 
the plant pole variation. Once the mappings have been 
plotted, it is not difficult to choose a minimum value 
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of P 1 such that the "far-off" closed loop poles are at 
the vertical boundary shown in Fig. 3.26 and a value 
of y such that the mapping of the plant pole varia- 
tion may be fitted inside the mapping of the dominant 
closed loop pole region for the same value of P 1 . 

The mapping of the dominant closed loop pole 
region with the mapping of the plant pole variation 
placed inside it is shown in Fig. 3.29. Since the 
mapping of the plant pole variation does not quite fit 
inside the mapping of the dominant closed loop pole 
region, a slightly larger value of y may have to be used. 

The values of kK, S and P are solved for 

o o 

using Eqs. 3.31. 3.32 and 3.33 and point A in Fig. 3.29. 
At point A 

X a = 800 Y a - 5450 

and from Fig. 3.26, at point A 



For y = 20000 and P. = 30, kK, S and P are 

1 o o 


kK 




= 800 - (- 6 ) ( 30) -10 = 970 


S 

o 


Y 

a 



1 


kK 


5450-10(30) 

970 


5.32 


FI0.9.Z9 CALCULATION Of SYSTEM CAIN AND 
COMPENSATION iE.RO POSITION FOR t » 20000 


Ht, 


IOOOO 


9500 


9000 


7500 


fcOOO 


5500 


5000 


+500 


4000 


3500 





8 ? 


p - — X- - Q . QP _ 20 6 i 

o kK 970 * 


Solving for the position of the compensation zeroes 
yields 

s = -2a 


a. - if - - 2 ’ 66 


n 2 2 

P„ = a +uj 
o z z 


uj = yp -a -a/20.61-5.32 
z o z 

- 3.60 


The position of the compensation zeroes Z, Z is then 


Z ,Z = -2.66 ± j3 . 60 

The actual closed loop poles for this choice 
of system gain and compensation zero position are 
found by determining the roots of the charar teristic 
equation 1+L^(s) = 0 where L^(s) is given by Eq. 1.12 
for points around the boundary of the plant pole 
variation. The convergence procedures in Appendix A 
were used to factor the resulting fourth order poly- 
nomial. The results of this operation are shown in 
Fig. 3.30 (irial 1). The points A and D corres- 
pond with those in Fig. 3.26. From Fig. 3.30, it is 
seen that the dominant closed loop poles lie outside 



Fib. 3.30 OOMfNANT ROOT T&ST FOR *« 20000 

AtfO /* ■ 30 

T'R.lA L JL • A K -- 9 70 £ * ~ Z . !»U +j* 3. 40 

TXfAL Z l kK-tG+5 j£» -Z.+U +j3.(*Z 
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their acceptable region for plant poles located at 
points B,B. The closest approach of the "far-off" 
closed loop poles also occurs when the plant poles are 
at points B,B which is -6.498 * J26.395- For all other 

points used in the dominant root test, the "far-off" 
closed loop poles were well t * the left of the 
established boundary shown in Fig. 3.26. 

In an attempt to place the dominant closed loop 
poles within their acceptable region, the position of 
the mapping of the plant pole variation within the 
mapping of the dominant closed loop pole region was 
slightly altered as shown in Fig. 3-29- For this new 
position and using the new coordinates for point A, 
the values of system gain and compensation zero 
position are 

k X = 1045 Z,Z = -2.46 ± j3 . 62 

The results of this trial are shown in Fig. 3.30 (Trial 
2). The dominant closed loop poles are still, outside 
their acceptable region when the plant poles lie at 
points B, B. The closest approach of the "far-off" 
closed loop poles is - 7.019 ± j 27 - 57 ^ when the plant 

poles ar~> at points B,B. Therefore the dominan closed 
loop pole region must be mapped into the X,Y pi me 
using a larger value of y. The mapping of the plant 
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pole variation for P 1 = 30 in the* £X, &Y plane 
remains unchanged. 

Figure 3*31 shows the mapping of the dominant 
closed loop pole region for y = 22000 with the mapping 
of the plant pole variation fitted within the interior. 
Using point A again gives the following 

kK = 1165 Z,Z = -2.432 i j3.608 

From Fig. 3.32, this value of gain and compensation 
zero position place the dominant closed loop poles within 
their acceptable region for all plant pole positions. 

The closest approach of the ”far-off” closed loop 
poles occurs again when the plant poles are at B,E 
and is -7.313 ± j 29 . 66 l . 

From inspection of Figs. 3.22 and 3.23. it is 
noted that the mapping of the plant pole variation is 
quite sensitive to the value of and that the size 
of the mapping of the plant pole variation decreases 
with decreasing P 1 . Therefore it may be possible to 
place the dominant closed loop poles within their 
acceptable region at a slightly smaller value of system 
gain kK if the value of P^ is decreased. This will 
also place the "far-off" closed loop poles closer in. 
Figure 3.33 slows the results of using a value cf P ] 

= 28 and the same of gain and zero position as for 
Trial 2 with v = 20000. The saving in gain between 


•> 1 


Flfr. 3.3/ CALCULATION OF SYSTEM GAIN AND 
COMPENSATION ZERO LOCATION FOR & * 22000 
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kK = 1045 and kK = 1 1 65 is approximately 1 db. The 

closest approach of the "far-off" closed loop poles 

is -6.324 ± j28.355 with the plant poles at B t B. 

This design is probably satisfactory since the 

closed loop poles lie to the left of the boundary 

shown in Fig. 3*26 for all but three other plant 

pole positions near point B used in the root test. 

From Figs. 3.29 and 3*31* it is noted that the 

mapping of the plant pole variation need not fit 

completely inside the mapping of the dominant closed 

loop pole region for all values of in order for the 

dominant closed loop poles to lie within their acceptable 

region. The requirement that must be satisfied is 

that the mapping of S s S in the AX, AY plane must 

J p i 

fit inside the mapping of the dominant closed loop 

pole region for S = S . From Fig. 3.31 f the Tine 
P p i 

segment AB corresponds to = -6 which is completely 
within the mapping of the dominant closed loop pole 
region for = -6. The line segment CD corresponds 
to S = 10 which is completely within the mapping of 

Jr 

the dominant closed loop pole region for = 10. 

So for in the design example, the problem of 
variation in the plant gain factor has not been 
considered. The actual value of gain which must be 


added to the system is given by Eq. 3.34, For the 




i 

i. 
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design of P 1 = 28, Z ,Z = -2.46 ± j3.62 and kK. = 1045, 

the value is 1045/k . . If the plant gain variation 

min 

is assumed to be 

k . = 1 £ k £ 1000 = k 

nun max 

the compensation networks must have a gain of at least 

1045. At k = k , the system gain is 1.045*10^ and 
max 

Tor this large value of gain, the dominant closed loop 
poles are essentially at the position of the compensa- 
tion zeroes. Using the method outlined in Section 10 
of this chapter, the angles of departure of the root 
locus for dominant closed loop poles located at points 
A, B, C and D in Fig. 3.33 were computed. The results 
were as follows: 

Point A: $ ^ = -14° 

B: = 69° 

C: <J d = 64° 

O: $ d = 133° 

These values are certainly satisfactory for the accept- 
able dominant closed loop pole region shown in Fig. 
3.33. Therefore the design is complete. 

3.12 Summary 

The design procedure for *_ain and plant pole 
variation is summarized below. 

1 . ) Map the acceptable dominant closed loop 
pole region into the X,Y plane using 


Eqs. 3.1^ and 3.15 with parameters y, 
and S • Map the plant pole varia- 
tion into the AX, AY plane using Eqs, 

3.16 and 3.17 with parameter . 

Determine the value of y and P^ such 
that the "far-off’ 1 closed loop poles 
just satisfy the minimum damping factor 
specifications for the problem and the 
mapping of the plant pole variation in the 
AX,AY plane can be fitted into the inter- 
ior of the mapping of the dominant closed 
loop pole region in the X,Y plane, 

2. ) Solve for the values of kX, S and P 

using Eqs, 3.11, 3.12 and 3.13 which 
determine the value of system gain and 
the c ompensat ion zero position. 

3. ) Check the design by determining the actual 

closed loop poles for plant poles lying on 
the boundary of the plant pole variation. 
This is accomplished by determining the 
roots of 1+L^(s) = 0. Check the final 
design for plant gain variation by using 
the method outlined m Section 10 of this 
chapter. 
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This design procedure which includes the effect 

of the "far-off" pole has several advantages over 

4 

the third order system considered by Horowitz. Using 
the third order approximation, the "far-off" closed loop 
poles can not be positioned so that the damping factor 
specifications are just satisfied. This results in 
a waste of system gain since the "far-off" open loop 
pole must then be placed sufficiently "far-off" so 
that its effect on the dominant closed loop poles is 
negligible. In the fourth order approximation, the 
effect of this "far-off" pole is considered on both 
the mapping of the dominant closed loop pole region and 
plant pole variation and can be used to determine a 
more economical design. 

The next chapter considers the additional 
effect of a drifting zero on the real axis on the 
design procedure. 








CHAPTER IV 


PROBLEM OF SIMULTANEOUS PLANT GAIN, 

POLE AND ZERO VARIATION 

4 . 1 Problem Definition 

This chapter presents an approximate design 
procedure for handling the added problem of a drifting 
zero on the real axis. This zero can be considered to 
be part of a plant with a transfer function P(s), 
given by 

p( s ) c - (4.1) 

s(s +S s+P ) 

P P 

where -z is the position of the drifting zero on the 

real axis. Typically this zero on the real axis is 

close to the origin and hence to the dominant closed 

loop poles, and has an appreciable ffect on the position 

of the dominant closed loop poles. The problem is 

then to choose the compensation zero position and the 

value of system gain such that the dominant closed 

loop poles lie within their acceptable region despite 

variations in the plant zero, parameter z, and in the 

plant poles, parameters S and P . The problem is 

P P 


shown in Fig. 4.1 
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4 . 2 Dosig it i J r 1 1 1 -r s o pi i y 

Any drifting zero of f(s? will appear as a 
drifting zero of T(s) (the system* t ranmi ssi ory, sirice 
i t la impossible to precisely cancel such zeroes. To 
minimize the effect of a drifting zero on T(s). a pole 
of L^(s) is placed near the zero in an. attempt t;> at 
least partially cancel the effect of this zero. If the 
system gain is sufficient ;Jv high, a closed loop pole 
will be \cry close to the zero despite its drift. The 
resulting dipole will have a negligible effect on the 
system response, if the system gain is large enough* 


4*3 T he Effect of Zero Variation on the L ami nan t 
Closed Loop Poles 


The effect of a poie^eru' pair on the real axis. 


near the dominant pole region, is shown in Figs. 4*2 
and 4,3. These figures are for a system designed 
according to the procedure in Chapter ill, with the 
celled pale-zero pair on the • real axis. The r.&ppiug:& 
vrere obtained by factoring the roots of Ub(?)=O t 
where L^is? is given by 


L .{ 3 

a ■ 


M 


2 .. - rj - w 

; + 3 o* + 1 V' S ; 


( Si *3 s*P ) t StP. 

P P 1 




and z and P ■ represent the poie-zeru pair on the real 
axis. The region of plant- pole variation is the same 
as was used in Chapter III, as is the acceptable dcmi- 
nant closed loop pole region. 
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Figure 4.2 illustrates the effect on the domi- 
nant closed loop poles for P =2 and z=3; it is similar 

z 

to the effect of lag compensation. Figure 4.3 is for 

P =2 and z=1 which has an effect similar to lead com- 
z 

pensation. These two figures, then, illustrate the 
effect on the dominant closed loop poles of a zero 
drifting from -1 to -3 when an open loop pole of L^(s) 
is placed at -2. As seen from Fig. 4.2, the system 
designed according to the procedure in Chapter III is 
not adequate to handle this zero variation on the real 
axis, since the dominant closed loop poles lie outside 
their acceptable region. The uominant closed loop poles 
could be forced into their acceptable region by using a 
larger value of system gain. However, this may result 

q 

in a waste of system gain. Horowitz presents a method 
to determine whether a design for a given region of 
plant pole variation is adequate to handle the effect 
of a drifting zero on the real axis for a specified di- 
pole (closed loop pole-zero) separation. The presence 
of the pole-zero pair on the real axis may influence 
the choice of the compensation zero position. The 
design procedure presented in this chapter takes into 
account the zero variation in determining the position 
of the compensation zeroes . 
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4 . 4 Design Equations 

In considering this problem, it does ret take 
long to determine that the drifting zero added to the 
plant pole variation, complicates the design procedure 
to a considerable extent. The dominant loop trans- 
mission must now be designed to take into account tvo 
independent types of variation, i,e., zero variation 
along the real axis and plant pole variation in the 
complex plane. Because of this, the nearest " far-off " 
open loop l'j is omitted from L^(s) to retain a fourth 
order representation fur the system. The expression for 
L d (s) is from Eq. 1,1 4 

kk(s 2 + S s+P ) (s+z) , kkn , (s) 

T / \ _ x o o' x ' is d v ' 

d s(s 2 +S a+P )(s+P ) d d^ 

v P P z 

(4.3) 

The expression for T d (s) is from Eq. 1.15 

P r p f p c / z (* +z ) p r P f P c A( s+Z ) 

r£ ( S } = 1 S — ^ Z 

d (s 2 +S r s+P r ) (s+p c )(s+p f ) 

Z 1 

(4.4) 

The real axis root locus for this system for the ex- 
treme positions of the drifting zero is shown in Fig. 

k.k. 

The mapping equations for this type of system 
are derived in the same manner as in Chapter III. The 
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characteristic equation of the system is, from Eq. 2.3 


d h( s ) = d H (' i ) +kKn ri( s ) 


(s 2 +S r s+P r )(s+p )(s+p f ) = s(s 2 +S s+P )(s+P z )- 

Z 1 F F 


kX( s +S q s+P o ) ( s+z ) 


or 


S +( s r +Pf + P C )s +(p f P c +S r (Pf + P C ) +p r )s + 

1 z 1 z 1 z 


<V p f/ p c )tS r p I p c ) S+P r p f p c 
1 z 1 z 1 z 


s^+VP +S +kh)s 3 +(P S +P +kX(S +z))s 2 +(p P +kh.(P +S z))s 
zp 7 x zpp v o v zp x oo /7 


+kKP z 
o 


(4.5) 


Equating the coefficients in Eq. 4.5 yields the follow- 
ing set of equations 


S +p_ +p = P +S +kK 
r 1, c z p 

1 z 


P f P c +S r^ P f . +p c ) +P r = P z S p +P P +kK(S » 
I z 1 z ^ 


+ z 


r (p„ +p )+S p p = P P +kh.(P +S z) 
f, c 7 r r f. *c z p ' o o' 

1 z 1 z v 


P p„ p = kKP z 
r*i . *c o 

1 z 


(4.6) 

(4.7) 

(4.8) 

(4.9) 
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For the problem considered in Chapter III, the 
parameters relating to the non-dominant closed loop 
poles on the left hand side of these quations were 
eliminated by substitution. In this case a slightly 
different approach is used. Define 


U = P +S +kX 
z p 

V = kXP z 

o 

X = P S +f +kK( S +z) 
z p p v o 7 

Y = P P +kx(p +S z) 

z p v o o 

The mapping equations for the 
region are then 


(4.10) 

(4.11) 

(4.12) 

(4.13) 

dominant closed loop "ole 


U = 


V = 


X = 


Y = 


S r + Pf 1 + P c 

1 z 


P r p f , p c 
1 z 

p f p c * S r (p f * p c >* P r 
1 z 1 z 

P r< p f/ p c >* S r p f p c 

1 Z I z 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


The total variation in X, AX, and The total variation in 
Y , AY, due to the variation in plant poles, i.e. para- 
meters 3 and P and the variation in the zero / are 
P P 

given by 





AX = P AS +AP +kXAz 
z p p 


(4.18) 


AY = P A P +kkS Az (4.19) 

z p a 

Since kh arid S q are not apriori known, the first 
approximations to aX and uY are taken to be 

AX = P AS +AP (4.20) 

z P p 

AY = P AP (4.21 ) 

z p 

These equations are now identical to the mapping equa- 
tions 3.16 and 3.17 in Chapter III. 

4 . 5 Design Procedure 

An outline of the design procedure that will be 
followed in this problem is as follows: 

1. ) Map the acceptable dominant closed loop 

pole region into the U,V plane using Eqs. 

4. 1 4 and 4.15 for fixed values of p and p . 

1 1 C Z 

2. ) Map the acceptable dominant closed loop 

pole region into the X,Y plane using Eqs. 

4.16 and 4.17 for the same fixed values of 
p and p c . 

1 Z 

3*) Map the plant pole variation into the 

uX,AY plane using Eqs. 4.20 and 4.21 for 

fixed values of P . 

z 

4.) Compare the two mappings in (2.3) above. 

If the plant pole variation in the &X,\Y 
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plane does xiA, fit within the interior of 
the mapping of the dominant closed loop 
pole region in the X,Y plane, then the 
dominant closed loop pole regi n will have 
to be mapped into both the U,V plane and 

X,Y plane using a larger value of p p . 

1 1 

5. ) Solve for the values of kh , S and P 

using Eqs. 4*10-4.13 where the values of 
U,V are obtained from the mapping of the 
dominant closed loop pole region in the U , V 
plane and the values of X and Y are obtained 
from the positioning of the mapping of the 
plant pole variation within the interior 
of the mapping of the dominant closed loop 
pole region in the X,Y plane. 

6. ) Determine the additional variation of the 

mapping of the plant pole variation in the 
X y Y plane by using Eqs. 4.18 and 4.19 and 
the values of kX and S obtained in ( 5 ). 

If this additional variation can not be 
accommodated within the interior of the 
mapping of the dominant closed loop pole 
region in the X,Y plane, then the dominant 
closed loop pole region must be mapped into 
both the U,V plane and X,Y plane using a 
larger value of . 


i 


I I o 


i 

i 
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f 
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The next five sections of this chapter elaborate on 
these steps in the design procedure, 

4 . 6 Mapping of the Dominant Closed Loop Pole Region 
Into the U,V Plane 

The design procedure begins by mapping the 

dominant closed loop pole region iuto the U,V plane 

with parameters p and p using Eq. 4 . 1 4 and 4.15- 

x ~ c 

1 z 

The parameter p„ , which is the ’'far-off” 

1 1 

closed loop pole on the real axis, plays the same part 

as the parameter y = kkP o did in the design procedure 

of Chapter III, i.e., large p implies large kk. 

1 

This can be ascertained from Eq. 4.6 since the para- 
meters S , p , P and S do not have a large variation, 
r * c z p ° 

z * 

The other parameter in this mapping operation 

is p , the closed loop pole near the drifting zero, 
z 

As an approximation to the value of p in the mapping 

z 

operation, it is assigned values that include the mini- 
mum and maximum vaiues of the zero z as well as inter- 
mediate values. In practice, this approximation can be 

improved by performing the mapping operations for p = 

'z 

Z max + ^1 and p c =Z min + 6 2 where ^ -j and - 9 are tde esti- 
z 

mated dipole separations when the zero is at its maximum 
and minimum positions respectively. The mapping of the 
dominant closed loop pole region into the U,V plane using 
Eqs. 4.1 4 and 4.15 is. shown in Fig. 4.5. 
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Fig. 

4.6. The 

unshaded 

region shown 

in Fig. 4.6, Ue 


noted by the set can be written as 


Q s iUi 


( 4 . 22 ) 


bow any point in will map into the acceptable dominant 

closed loop pole region in the s -plane for all values of 

o = p,. C r z . ,z 1 . This is the same tvpe of ^ r^ru- 
c -zi min max- ' ^ 

merit used in Section 7 of Chapter III. An arbitrary 

point in w will he used later to solve for the required 

val u*> s o £ s vs t e rn gain kh a nd S , P wn ion d e t e r ahne 

o u 

the compensation zero position. 

1 n order to insure that rh o exists, i.e., 

Q / * where is the null set, restrictions LtiU^ i be 
placed or. the magnitude of the zero variation on the 
real axis and hence p 

c 

z 

consider the mapping of a point in the dominant 


closed loop pole region characterized 
for p -p ♦ Equations 4.1 4 a: d 4 

Z zmir: 

into the U,Y plane for a fixed value 


dv ana IV, 

1 1 - 1 
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of p ^ as 1 O 1 1 1 ; V S l 
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z zmax 


same fixed value of p« , i.e. 

■1 


* r +P f +? c 
2 1 zmax 


p f Pc. 


= V, 


(4.25) 

(4.26) 


zzaax 2 

If the mapping defined by Eqs. 4.23, 4.24 and Eqs. 
4.25, 4.26 are to have at le&st one point in common, 
the following conditions must be satisfied; 



= u 2 ; v 

1 2 V 2 

• S 

r 2 

* p c ~ p c 

zmax zmin 


' Z. Olii X 
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znxn 
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order to 
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(4.28) 


th^re must exist two points contained in the acceptable 


dominant rinsed loop pole region defined by S , P 
» r 1 5 


xnd £> , P„ respectively such that Lqs. 4.27 u?.« 


4.28 are satisfied. These equations set a limi; on the 


allowable variation' in p and hence z in t ti‘j 

c 

z 

mapping cp.tition. Thssv restrictions ere not overly 


stringent as will be illustrated in the design example 
in this chapter. 
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If the zero variation is large compared to the 
acceptable dominant closed loop pole region and Eqs. 

4.27 and 4.28 are not satisfied, the design procedure 
presented in this chapter *ould have to be modified to 
obtain a design. 

The effect of the parameter p for fixed p^ 

c z 1 1 

on the mapping of the dominant closed loop pole region 
into the U,V plane is shown in Fig. 4.7- The acceptable 
dominant closed loop pole region used in this and sub- 
sequent mappings in this chapter is the same as that 
used in the design examples in Chapters II and III and 
is shown in Fig. 3*26 and in the root tests of Chapter 
III. 

The effect of the parameter p r at fixed p 

1 c z 

on the mapping of the dominant closed loop pole 
region into the U,V plane is shown in Fig. 4.8. 

4 . 7 Mapping of the Dominant Closed Loop Pole Region 
Into X.Y Plane 

Th6 mapping of the dominant closed loop pole 
region into the X,Y plane is accomplished using Eqs. 

4.16 and 4.17 with the sane values of p- arid p 

X 4 C 

1 z 

that were used in the mapping in the previous section. 

This mapping operation is shown in Fig. 4.9. As in 

Chapter III, the mapping of the plant pole variation 

will be fitted inside this mapping to solve for the 

values of kk, S and P . 
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The effect of the parameter p for fixed 

z 1 

on the mapping of the dominant closed loop pole region 
fnto the X,Y plane is shown in Fig. 4 C 10. The effect 
of the parameter p for fixed p on this same mapping 

I 4 C 

1 z 

is shown in Fig. 4.11. 

4 .8 Mapping of the Plant Pole Variation Into the 
AX, AY Plane 

The napping of the plant pole variation into 

the AX, AY p? > . ^ accomplished using Eqs. 4.20 and 

4.21. Since equations are identical in form as 

those in Chapter III, they will not be investigated in 

detail here. This mapping operation only considers 

the variation in the plant poles. The variation in the 

zero on the real axis will be taken into account later 

in the design. The variation in plant gain may be 

handled in the same manner as outlined in Chapter III. 

The parameter P r , the open loop pole placed 
z 

near the drifting zero, occurs in the same manner in 

Eqs. 4.20 and 4,21 as P did in the mapping of the plant 

pole variation in Chapter III. 

The mapping of the plant pole variation into 

the AX, AY plane using Eqs. 4.20 and 4,21 is shown in 

F^‘g. 4.12 for various values of P r . The plant polo 

z 

variation used in this mapping is identical to that 
used in Chapter III and is shown in Fig, 3.26. As in 
Chapter III, the units on this mapping in the AX, AY 
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plane must be the same as those used in the mapping of 

the dominant dosed loop pole region in the X,Y plane. 

The mappings of the plant pole variation for P r =2 and 

z 

comparable units as those used in the mapping of the 
dominant closed loop pole region into the X,Y plane 
are also shown in Fig. 4.12, 

4 . 9 Calculation of System Gain and Compensation Zero 
Location 

Once the mappings of the dominant closed loop 
pole region are obtained, an attempt is made to fit 
the mapping of the plant pole variation in the AX, AY 
plane into the interior of the mapping of the dominant 
closed loop pole region in the X,Y plane as shown in 
Fig. 4.13* If this is not possible, then the mapping 
of the dominant closed loop pole region must be per- 
formed at higher /ajues of p^ . The mapping of the 
plant pole variation in the AX, AY plane is unchanged 
since is assumed to be fixed at the beginning of the 

design. Once the mapping of the plant pole variation 

- ♦ - * 

can be fitted inside the mapping of the dominant closed 
loop pole region in the X , \ r plane, the value of system 
gain kk and S , P , trie parameters that determine the 
compensation zero position, can b< determined. 

The first step in solving for kk , S anc P is 
to choose an arbitrary point within the set Q of the 
mapping of the dominant closed loop pole region in the 
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U , V plane* Denote this point by "0 U and the values 

of U and V at this point as U and V . 

o o 

At point A in Fig. 4.13, denote the values of 

S P , X and Y as S- * P- * X and Y^ respectively, 
p p Pa P a a a 

Since P is also known, kX can be found from Eq. 4.10 
z 

i.e . , 

kh = li -P -S (4.29) 

ozp 

This leaves Eqs . 4.11, 4. 12 and 4.13 to solve for 

and 5 , The value of a is also an unknown. This value 
o 

of 2 is needed to compute the added variation in the 
mapping of the plant pole variation in the X,7 plane 
(Eqs. 4.18, 4.19)* The effect of the zero variation 
on the -design is considered in the next section. 

Using Eq. 4.11 to eliminate z in Eq. 4.12 and 
4.13 yields 

X = P S +P +kK(s +V /khP ) (4.30) 

a z p d: ' o o' o ' ' 

*a *a 

Y„ = P P +kk(p +S V /k_kP ) (4.31) 

a z p 'ooo o - ' 7 


Solving for in these two equations yields 


ing third order equation for 

(P P -Y )P 2 V (X -P S -P 
o3 2 P a a ° ° a Z P a P a 


)P 


kk 


(kk)‘ 


(kK)‘ 


the follow- 



(4.32) 



The design xs not complete once these values 
have been obtained since the zero variation Az in 
£qs. 4,18 and 4.19 was neglected in the mapping of the 
plant pole variation into the AX, AY plane. This added 
variation in the mapping of the plant polo variation is 
taken into account in the X,Y plane by considering that 
the coordinates of any point on the mapping of the 
plant pole variation in the X,Y plane could change by 
as much as 


AX = kKAz 
z 


(4.34) 


AY = kKS Az 
z o 


(4.35) 


where AX and AY is the change in any point on the 
mapping of tie plant pole variation -in the X,Y plane 
due to the variation in the drifting zero only. The 
zero variation Az is computed as follows: From Eq. 


4.1 I, the nominal value of z denoted by z is 






(*.36) 


12 



o 

This is the value of z when the mapping of the plant 

pole variation in the X,Y plane is in its original 

position. (The position used in the previous section to 

compute k 5 and P ). The nominal value of z, z t 
o o' o 

will lie somewhere in the interval fz’ ,z ] because 
of the values of p c ^ used in the mapping operations. 

The positive change in the coordinates of any point 
on the mapping of the plant pole variation in the X,Y 
plane is 


AX + = 4cK(z -z ) = 
z x max o ' 

A Y + = kXS (z — z ) = +kKS A z 

z o' max o' o 

whereas the negative change is 


(*.37) 

(*•38) 


AX Z = kK( Znjin -z o ) = -kXAz (*.39) 

AY~ = kXS (z . -z ) = -kKS A z (*.*0) 

* o' mm o' o v ' 

This variation is shown in Fig. *.13. Note that AX + 

z 

is not necessarily equal to AX~ since z may lie 

z o 

anywhere in th«- interval fz . ,z ]. This added 
variation is taken into account by considering the 
changes in the points A, B, C arid D on the boundary 
of the mapping of the plant pole variation in the X,Y 
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plane as shown in Fig. 4.13. The maximum possible 
plant variation is the figure A'A'B'C'C'D* where the 
primes indicate the new positions of the points A, B, 
etc. If this added variation can not be fitted within 
the interior of the mapping of the dominant closed loop 
pole region in the X,Y plane, then, the dominant closed 
loop pole region must be mapped into the U,V and X,Y 
planes using a larger value of Pf ^ . The mapping of the 
plant pole variation into the AX, AY plane remains un- 
changed as long as P z is unchanged. 

The design should be checked when the total 
Variation can nearly be fitted within the mapping of 
the dominant closed loop pole region in the X,Y plane 
since this is an approximate design procedure. The 
assumption that p^ remains constant as the dominant 
closed loop poles vary is only approximately true. 

The mapping of the dominant closed loop poles at fixed 
Pf.j essentially determines the minimum value that p^^ 
attains when the actual closed loop poles are obtained 
by factoring l+L^s) = 0. The mapping of the dominant 
closed loop pole region in Chapter III at constant 
Y = kKP Q was completely valid since kK and P q are 
fixed if the gain variation is neglected. 


4.11 Summary 

Unlike the design procedure developed in Chapter 
III, the design procedure developed here is approximate. 
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The assumption that remains fixed as the dominant 

closed loop poles vary is not strictly valid* For 
systems with large plant parameter variations, this 
approximation should lead to an acceptable design since 
large gain Implies a large value of p^ and thus its 
effect on the dominant closed loop poles wixl be slight. 

The major difficulty encountered in this problem 
was determining what effect the zero variation had on 
the mapping of the plant pole variation (Eqs. 4.18 and 
4.19)* From these equations, the values of kX and 
were needed to map the plant pole variation into the 
AX, AY plane exactly but this mapping itself was needed 
to solve for kK and S^. This difficulty was overcome 
by first considering the effect of the plant pole 
variation alone and then checking to determine if this 
design was adequate to handle the added zero variation. 
The effect of the zero variation on the mapping of the 
dominant closed loop pole region was taken into-account 
by the parameter p 

c z 

4.12 Design Example 

The design example presented here has the same 

3-domain specifications as the design example iri Chapter 

III (Fig. 3*26) with these two exceptions: The added 

zero variation along the real axis, i.e. z ^z-z 

mm max 

where z , =1 and z ^3 * The effect of the "far-off" 
min msx 
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pole P.| is neglected in order to retain a fourth order 
representation for the system. 

The first step in the design is to choose a 
fixed value of . In this design procedure, it is 
difficult to obtain art approximation for the value of 
Pfj that should be used. Obviously is going to be 

relatively far-removed from the acceptable dominant 
closed loop pole region for large parameter variations. 
Also the time domain specifications for the problem 
will probably dictate some minimum valvie of pf^ (See 
Fig. 4.1). A value of 40 for p^ was chosen as the 
first estimate. Figures 4.14 and 4.15 illustrate the 
mapping of the dominant closed loop pole region in the 
U,V and X,Y plane respectively for p^ = 40. In 
addition, the dipole separations for z , 6., and z . , 

6 2 ,have been estimated at 0.3 and 0.05 respectively 
(See Section 4.6). Only the mappings for the maximum 
and minimum values of p Cz are shown in Fig. 4.1 4 
since the set Q is completely defined by these two 
mappings . 


A value of P =2 will be used since this is mid- 
z 

way between the extreme zero positions and is tho value 

used to illusr.rate the effect of zero variation on the 

dominant closed loop poles (Figs. 4.2 and 4.3). The 

mapping of tho plant pole variation for P =2 into the 

z 

AX. AY plane is available in Fig. 4.12. 
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I'i.guro 4,1 5 shows the mapping »f th«; pla it 

pole variation for P =2 within the mapping of the 

z 

dominant closed loop pole region for =40 and Pc z 
= 0.95, 2 and 3.3. At point, A in Fig. -4,15 


X * 420 

a 


= 2025 

a 


From Fig. 3.26 at point A 


S ^ -6 

P a 


F =10 
Pu 


From Fig. 4,1*4, point "O’* in the set b is arbitrarily 
chosen t» have the coordinates 


U = 31 

o 


V =2100 
o 


Using P =2 , the value of kK can be obtained from Fq« 
z 

4.29, i-e. 


kK = U -P -S = 51-2+6 
o z 


= 55 


Using Fcj , 4.32, the third order equation for \* ^ is 

(P P -Y ) V U -P H -P )p 2 

« Z P, a , o a Z p p u V 
_ 3 •*- .+2 4 a 1 a o 

P + -rr F + - = o 

o Fh o / , , + 2 z , , v 2 


(kK)' 


(kkf 


,3 + (2(lO)-2Q ^) r 2 jjMOOj; _ 

T »• J r . — 


>2 


o 


(35) 


( 55 


J ) 


r j -'i6. 5rV 295P - 1 460 =: 

n ' ^ t) O 


0 


o 


13 ^ 


Factoring the one real positive root of interest using 
the convergence procedures in Appendix A yields- 


P q = 27.85 


is obtained from Eq. 4.33, i.e. 


X--P S n -P D V 

a z Pa Pa o 

kX ~ kXP 


420 - 2 (- 6) -10 
55 

= 6.-31 


■ 2100 

55 ( 27.857 


The position of-the compensation zeroes is 

=-^= - 3.105 


• /D 2 2 

*2 = ^ P o' a 2 


= W18.20 = x4 .266 


The design can now be checked to determine if it is 
adequate for the zero variation* The nominal value 
of z , z , ^5 Tound from Eq. 4.36, i.e. 


z = 
c 


V C 2100 _ . 

kKP Q “ 55(27.85 ) “ '* w 


The positive variation in the mapping of the pl^nt pole 
variation in the X, T plane is found from Eqs. 4.3/ 
and 4.3S. i.e. 


eX* 5 = tdv{z — z } = 55 ( 3 . 00-1 . 37 ) 
z ' max o r * 


« *89.? 
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AY + = kkS (z -z ) 
z o' max o 


= 55(6.30(1.63) = +566 


The negative variation is found from Eqs. 4.39 and 
4 . 40 1 , i . e . ' . 


AX = kk (z . -z ) = 55(1 .00-1 .37) 

z ' min o' 


= - 20.4 


AY~ = RKS (z . -z ) 
z o' mm o' 


= 55(6.31 ) (-0.37) = -128.5 


This added variation is taken into account by consider- 
ing the change in the coordinates of points A, B, C and 
D in the X.Y plane as shown in Fig. 4.15. Since this 
added variation can not be completely accommodated 
within the mapping of the dominant closed loop pole 
region, the design is probably not adequate. Using 

these values of kK, S and P , the actual closed loop 

o o 

poles were found for the boundary of the plant pole 
variation shown in Fig. 3.26 by determining the roots 
of 1+L^(s)=0. The results are shown in Fig. 4*1 6 
which confirms that the design is inadequate. 

As a check on the design procedure . the posi- 
tion of the compensation zeroes was arbitrarily moved 
t o 


z,z = -3-75 ± j4 . 25 


FlO. +.11* COM INANT ROOT TEST FOR 7 %*+0 

ANO /J ■ t 
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in an attempt to force the dominant closed loop poles 
into their acceptable region. The result using the 
same value of gain as before is shewn in Fig. 4,17 
which indicates that this compensation zero position 
does not ^yield a satisfactory design either. This 
implies that the mappings must be performed using a 
larger value of . 

The next value chosen for p r was 50. The 

*1 

mapping of the dominant closed loop pole region into the 

U,V plane aid X.Y plane for this value of p~ is shown 

in Figs. 4.18 anc t 4.1^ respectively. The same value of 

P =2 was used so the mapping of the plant pole varia- 
z 

tion in the AX, AY plane is unchanged. The mapping of 
the plant pole variation fitted within the interior of 
the mapping of the dominant closed loop pole region in 
the X,Y plane is also shown in Fig. 4.19. At point A 
in Fig. 4.19 


X a = 525 Y a = 2500 



The point "O'* in the set Q from Fig. 4.18 is arbitrarily 
chosen as 

U o = 61 V = 2650 

o 

Using P =2 , the value of \cK from £q. 4.29 is 65 . Using 
z 

Fq. 4.32, the third order equation for p is 
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P 3 - 38.18 P 2 + 331 P - 1663 = 0 
o ' o o 

The real root of interest is 


P =28.53 

o 

S from Eq. 4. 33 is 
o 

S = 6.68 
o 

The position of the compensation zeroes is 

S - -3.3<* 

= ± I >.17 

The nominal value of z, z q , from Eq. 4.36 is 

z =1.43 

o 

The added variation in the mapping of the plant pole 
variation is found using Eqs. 4.37 - 4.40, i.e. 


AX + = +102 
z 


AX~ = -?8 
z 


AY + = +682 
z 


AY - s -186.5 
z 


This additional variation is also shown in Fig. 4.19* 
The additional variation can be easily accommodated 
within the mapping of the dominant closed loop p »le 
region. The actual closed loop poles for this value 


of system gain and compensation zero position arc shown 


1 42 

in Fig. 4.20* The design is more than adequate to handle 
the plant pole and zero variation. The difference be- 
tween the gain of this design and the previous design 
is approximately 1.7 db. 

For this particular design, a system gain of 
less than 65 would probably be adequate. This could 

be verified by using a value of p^ = 45 for the design 

*1 

procedure. This would result in a system gain of 60. 

Further reduction in system gain may be possible by 

varying the value of P , the open loop pole used to 

z 

partially cancel the effect of the drifting zero. 

As can be anticipated from the root test shown 
in Fig. 4.20, the variation in plant gain will not 
result in the dominant closed loop poles leaving their 
acceptable region. Using the procedure presented in 
Chapter III to check the angle of departure of the 
dominant closed loop poles for the root test shown 
in Fig. 4.20 confirms this. 



F/fr. 4.20 DOMINANT HOOT TEST FOR • SO 

AND ^-2 

ktC'OS £• ~ 3.3+ +j +./7 


DOMINANT CLOSED 
LOOP POLE REGION 






CHAPTER V 
CONCLUSIONS 

The design procedures presented in this paper 
are for fourth order systems with large variations in 
the plant parameters (gain factor, poles and zeroes) . 
Extension of these design procedures to a fifth order 
system would be difficult since the mapping equations 
would be very cumbersome. 

Two possibilities for additional work in this 
area are presented in this chapter which would be a 
significant improvement over the design procedt*res 
developed in this paper. 

The first of these is a computor design routine. 
An initial guess would be made as to the compensation 
zero location at a fixed value of system gain.. By 
examining the resulting dominant closed loop poles and 
possibly using a gradient technique, one could determine 
in what direction the compensation zeroes should be 
moved to place the dominant closed loop poles within 
their acceptable region. This is complicated by the 
fact that the necessary value of system gain is also 
unKnovn since for small values of system gain, a design 
is impossible for any compensation zero location. 



1 U J 

Analytic expressions for the boundary of the plant pole 
variation and the acceptable dominant closed loop pole 
region would probably be required* 

The second possibility in this problem, would 
be an attempt to obtain an analytic solution for the 
compensation zero location and necessary value of 
system gain to place the dominant closed loop poles 
within a given region in the s-pluue for a given region 
of plant pole variation. Intuitively, one wauid think 
that there is a unique value of system gain and compen- 
sation zero location such that the closed loop poles 
lie within their acceptable region and the system gain 
is minimized, Unfortunately, a solution in closed 
form appears to be very difficult to obtain in even 
trivial cases. Analytic expressions for the boundary 
of the plant pole variation and the acceptable dominant 
closed loop pole region would certainly be required in 
this de s igri procedure. 
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APPENDIX A 
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POLYNOMIAL FACTORING 

A . 1 Statement, of Problem 

In many investigations of- feedback control 
systems, polynomials of rather high degree must be 
factored* In the analysis of feedback control sys- 
tems, the approximate location of the roots in the 
complex plane are often known. The problem usually is 
one of determining if the roots lie within an acceptable 
region despite variations in gain and/or plant para- 
meters. The polynomial considered in this appendix 
is of fifth order. A procedure is first developed for 
extracting one real root from this polynomial. The 
resulting fourth order polynomial is then factored 
into the product of two quadratics using Lin's 
method J This method is applicable whether the fourth 
order polynomial has real or complex roots. A method 
" is then presented for extracting two real roots froir. a 
fourth order polynomial. The nomenclature is r-hosen 
to aid in the programming of these procedures on a 
digital comput or. The convergence of any of these 
methods is not guaranteed. 
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A. 2 Extraction of One Real Root From A Fifth Order 
Polynomial 

The fifth order polynomial ic 

s 5 +¥1s i *+W2s 3 +W3s 2 +V4s+W5 (A.1) 

The coefficients ¥1 , ¥2 , etc., are considered to be 
real in all cases. 

This procedure is based on the fact that an 
approximation to the real root is given by the quotient 
¥ 5 /¥ 4 . An improvement on this approximation can be 
made after one long division trial by considering the 
binomial term in the last subtraction process. 1 ¥ r ith 
this in mind define 

v(j) = V5/\ik J^I only (a. 2) 

where: V(j) — J th approximation to the real root. The 

long division operation is shown below. 

s 4 +A1 ( j)3- ? +A2(j)s 2 +A3(j)s+A*t(j) 
s+V ( J ) j s 5 +Wl s +W2 s- 5 +W3s 2 +W4s+¥5 

s^-t-V ( J ^ s 4 
{ ¥ T - V ( J ) } s^+¥2s^ 

' fWI-V(j)] V(j)-V 2 ( J)] s'* 

{ ¥2 -W 1 V ( J ) + V 2 ( J ) } s 3 +W3s 2 

f ¥2-¥1 V(J)+V 2 { J)} s 3 +f¥2V( J)-¥1 V 2 ( J )+V 3 ( j)) s 2 
{W3-W2V( J)+¥1V 2 (J)-V- 3 (J)J s 2 +W4s 


{ W5-V4V ( J ) +W3V ( J ) -V2V J ( J ) +Vi 1 v ( J ) -V 0 ( J ) } 

The remainder term is the test for convergence. Le 

X = W5-W4v(j)+W3V 2 (j)-V2V 3 (j)+W1V 4 (vl)-V 5 (j) (A. 3) 

If X is not sufficiently small, then the next trial 
divisor V(j+l) is given by 

W4-V3V ( J ) +¥2 V 2 (j)-V1V 3 (j)+V 4 (j) 

The coefficients of the fourth order polynomial 


are given by 

A 1 ( J ) = W 1 —V ( J ) (£.5) 

A2 ( J ) = W2-VT 1 V( J) +V 2 ( J ) (A. 6) 

A3(«») = W3-W2V(J)+V1V 2 (J)-V 3 (J) (A.7) 

a 4 ( J ) = K4-W3V(j) + W2V 2 (j)-W1V 3 (j)+V 2f ( J) (a. 8) 


The iteration procedure is continued until the 
value of X is sufficiently small or until a sufficient 


number of trials have been made. 
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A. 3 Lir^s Method 

This method is essentially the same as the pre- 
vious procedure except the trial divisor is a quadra- 
tic term. Let the fourth order polynomial be given by 

s Zf +A1s 3 +A2s 2 +A3s+A4 (A ,'j) 

Define 

B1(I)=A3/A2; B2(i)=a 4/A2 1=1 only (A. 10a, 10b) 

where: B1 (i) , B2(l) = coefficients of I 13 * trial quadratic 

The trial divisor in the long division operation is 
given by 

s 2 +B1 (l)s+B2(l) (A. 11) 


Performing the long division operation results in the 
following 


s 2 +B1 ( I ) s+B2 ( I ) 


s 2 +C1 ( I ) s+C2 

T a o 

s +A1s +A2s~ 


ill 

+a3s+a4 


s*+B1 (l)s 3 +B2(l)s 2 

{A1-B1 (I )} s 3 +{A2-B2(l)} s 2 +A3s 


{ A1 -BltflkM A1 BlCO-BI 2 (l)fs 2 +{ A1 B2(l)-Bl(l)B2 (l)|s 
{ A2-B2 ( I ) — A1 B1 (l)+B1 2 (l)} s 2 + { A3-A1 B2 ( I ) 

+B1 ( I ) B2 ( I ) } s+A4 

{ A2-B2 ( I ) -A1 B1 (l)+B1 2 (l)] s 2 + [A2B1(J.)- 
B1 {I)B2(I)-A1B1 2 (I)+B1 3 (l)}s 
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+f A2B2 (i)-B 2 2 (I)-A1B1 (l) B2 ( I ) f B2 ( I ) B1 2 ( I ) } 

{ A3-A1 B2 ( I) + 2B1 ( I ) B2 ( I ) -A2B1 ( I ) +A1 B1 2 ( I ) 

-B1 3 (l)} s+{a4-A2B2(I)+B2 2 (i)+A 1B1 (l)B2(l) 

-r2(l)B1 2 (l)} 

The two remainder terms are the test for convergence. 

Let 

X=A3-A1 B2 (I) +2B1 (i )B2 (i) -A2B1 (i ) +A1 B1 2 (l ) -B1 3 ( I ) 

(A. 12) 

Y = a4-A2B2(i)+B2 2 (i)+A1B1 (I)B2(I)-B2(I)B1 2 (I) 

(A. 13) 

If X and Y are not sufficiently small, then the co- 
efficients of the next trial quadratic divisor are given 
by 

B1 (l+ 1 ) = A3-A^ll) ± Bl(llB2(|l_. U . l4 ) 

A2-B2 (i ) -A1 B1 (l)+B1 2 (l) 

B2(I+1) = — (A. 15) 

A2-B2 ( I ) — A1 B1 (l)+B1 2 (l) 

The coefficients of the remaining quadratic are 

Cl (l) = A1 -B1 (I) (A. 16) 

C2(I) = A2-B2(l)-AiB1 (I)+B1 2 (l) 


(A. 17) 
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The iteration procedure is continued until the 
values of both X and Y are sufficiently small or a 
sufficient number of trials have been performed. 

A. ^ Extraction of Two Real Roots From a Fourth Order 
Polynomia l 

This method is applicable when it is known that 
the fourth order polynomial has at least two real roots 
The fourth order polynomial is given by 

s^+A1 s 3 +A2s 2 +A3s+A4 

As with the fifth order polynomial, define 

V(J) = A4/A3 J = 1 only (A.18) 

where; V(j) = 3 th approximation to the real root 


The long division operation results in 


s+V ( J ) 


s 3 + Pl ( J)s 2 +D2( J)s+D3(J) 
4 q o 
s +A1s~+A2s +A3 s+A4 


sVv( j)s 3 
{ A1 ~V( j)} s 3 +A2s 2 


{ AT -V ( J ) ] s 3 +{A1 V( J)-V 2 ( J)} s 2 
{ A2-A1 V( J)+V 2 ( J)} s 2 +A3s 


{A2-A1V( J)+V 2 (j)} s 2 + { A2 V ( J ) - A1 V 2 ( J ) +V 3 ( J ) } s 
{ A3 -A2V ( J ) + A 1 v 2 ( J ) -V 3 ( J ) } s +A4 

{ A3-A2V(j)+A1 V 2 ( J)-V 3 ( J)} s+ { A3V ( J ) -A2V 2 ( j) 

- ■ + A 1v. 3 (J ) ,-v\ -U I 

{a 4-A3V( j)+A2V 2 ( j)-aiv 3 ( j) +v*( j)} 
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The remainder term is again the test for convergence, 
i.e. , let 

X = a4-A3V(J)+A2V 2 (j)-A1V 3 ( J)+V** (j) (a. 19 ) 

If X is not sufficiently s all, the next trial 
divisor is 

V(J + 1) = — — 5 = (A. 20 ) 

A3-A2V(J)+A1V*(J)-V J (J) 

The coefficients of the third order polynomial are 
D1 (J) = A1 -V(j) (A. 21 ) 

D2 ( J ) = A2-A1V( J)+V 2 (J) (A. 22) 

i>3(j) = A 3 -A 2 v(j)+aiv 2 (j)-v 3 (j) (a. 23 ) 

After the first real root is obtained with a sufficient 
degree of accuracy, the second real root is ex\ -icicd 
from the remaining third degree equation. Let the 
third order polynomial be defined by 

s 3 + C1 s 2 +C2s + C3 (A. 24) 


Define 


u(l) = C3/C2 1=1 only (A. 25) 

where: U(l) = 1 th approximation to real root. 

The long division operation results in 


1 r ,k 


s 2 +B1 s+B2 

s+U( I ) | s 3 + C1 s 2 +C2s+C3 

s 3 +U(l)s 2 

[Cl -U(l)} s 2 +C2s 

f C 1 -U ( I ) 1 s 2 + {C1U(l)-U 2 (l)l s 
{C2-C1U(l) +U 2 (e)} s+C3 

f C2-C 1 U ( I ) +U 2 ( I ) } s+{C2U(l)-C1U 2 (l)+U 3 (l)} 
f C3-C2U(l)+C1U 2 (l)-U 3 (l)} 

The test for convergence is 

X = C3-C2U(I)+C1U 2 (I)-U 3 (I) (a. 26) 

The next trial divisor, if X is not sufficiently small, 
is 

U(I+1) = — 5 (a. 27) 

C2-ClU(l)+U (i) 

The coefficients of the quadratic equation are 

B1 (I) = Cl -U(l) (A. 28) 

B2 ( I ) = C2-C1U(I)+U 2 (I) (A. 29) 

In some cases convergence fails in the case of extracting 
one real root from a third order polynomial. In this 
case convergence can sometimes be obtained by extracting 


a quadratic from the third order polynomial. 

Define the third order polynomial as 


s +C 1 s +C2S+C3 (A. 30) 

Also define 

D1 (I) = C2/C1 ; D2 ( I ) = C3/C1 1=1 only (A. 31a, 

A. 31b) 

where: D1 (l) , D2(l) = coefficients of the trial 

quadratic . 

The trial divisor in the long division operation is 


given by 


s +D1 (l)s+D2(l) 


(A. 32) 


The long division operation results in the following 

s+V ( I 


D1 (i ) s+D2 ( i) | s' , +C1s +C2s+C3 


s 3 +D1 (l)s 2 +D2(l)s 


(C1-D1 (I)J s +{C2-D2(T)} s+C3 

{C1-D1 (i)} s 2 +{C1D1 (l)-D1 2 (l)} 
+ { C1D1 ( I ) -D1 ( I ) D2 ( I ) } 


{ C2-D2 (l)-C1D1 (l)+Dl(l)J s 
+ {C3-C11)1 (i ) +D1 (l)D2(l)} 

The two remainder terms are a test for convergence. Le 


X = C2-D2 ( i) -Cl D1 (i)+D1 2 (i) 


fc.33) 


Y = C3-C1D1 (l)+D1 (I)l)2(l) 


(A.3M 


If X and Y are not sufficiently small, then 
the coefficients of the next trial quadratic are given 





1 5 0 


J)1(l + 1) = 

1)2 ( I = 1 ) = 

The real root is 
s = -V ( I ) 

A . 5 C onclusion 

The factoring procedures presented in this 
appendix may be used for fifth, fourth or third 
order polynomials possessing either real or complex 
roots or both. It should be stated again that the 
convergence of none of these procedures is guaranteed. 


C2-I)2(l 


Cl -D1 1 


C3 

C1-D1 (I) 


= - (Cl ~D1 (i) ) 


(A. 35) 


(A. 3 6) 


(A. 37) 


APPENDIX B 

ANGLE CONTRIBUTION THEOREM 

B.1 Statement of Problem 

This appendix presents a geometric proof of 
the statement in Chapter II regarding the angle contri 
bution of two zeroes to a complex pole located in the 
s-plane. The statement to be proved is as follows: 

The angle contribution due to two conjugate zeroes to 
a complex pole is a constant if the zeroes are located 
on a circular arc drawn through the complex pole and 
its conjugate and a third point X on the real axis 
defined by the equation 

/.xp* = 6 z /2 (B.1) 

where: p* is the complex pole 

6 z is the angle contribution due to the zeroes 

B . 2 Geometric Proof 

The geometric proof essentially shows that if 
the compensation zeroes are located on a circular arc 
through the points * the angle contribution is 

a constant independent of the position of the zeroes 
on the circular arc and that Eq. B.1 is satisfied. 

The proof is shown in Fig. B.1. 



The proof begins by defining 


(B.2) 


*1 - L z ps 


(B.3) 


*2 - 1 2 Pd 


to prove that 


4' 1 1-4* 0 = 2a 


( B . 5 ) 


From Fig. B.l, it is noted that the inscribed angles 
p*Zp*, p*Xp* f and p*Zp* all intercept the same circular 
arc p*p*. Therefore from a basic theorem in plane 
geometry, 1 the following can be stated 


L Pd Z Pd = V2m p d p- 


(B.6) 


L Pd X Pd = V < Zm Pd p d 


(E.7) 


L PgZPS - ¥2m p d p* 


(B.B) 


where: 1/2 tn o,p*; = one— half ox'' the measure of th< 

* d*a 

circular arc p*D* 

Also from Fig. B.1 f 


/. p*Zp* = <*> 

a d 1 c 


(B.9) 


L p*Xp* = 2a 
*a u 


(B.10) 


£ Pd Zp d = 


(B.1?) 


Therefore 


^ j -^2 = 2a 




